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Silo is a common storage structure for storing 
granular materials like food grains. Reinforced concrete 
silos are preferred over other convensional storage stru- 
ctures because of various reasons, e.g., economy, saving 
and better preservation of grains. 
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The design for a R.C. Silo is essentially a 
labourious trial and error process and cannot be repeated 
many times by hand computation, thus resulting in general 
m an uneconomical structure. The forces acting on a silo, 
e.g. , the dead load, the wind load, the earthquake force, and 
the grain pressures depend directly on the shape and dimen- 
sions of the structure. Thus the attempt of an economical 
design and saving m constructional materials couples itself 
with the reduction of these forces. Any saving achieved m 
the cost of R.C. Silo is going to build-up because of the 
repetitive requirements of this class of structure. Hence, 
the problem is most suitable as an automated optimum design 
problem. 

A cast-in-situ circular R.C. Silo with conical 
bottom, supported on a system of rmg-girder and columns 
has been considered m the present work. A mathematical 
programming problem has been formulated on the basis of 
Indian Standard Specifications with minimum cost of the 
material used as the objective function. The SUMT algorithm 
has been used to seek the optimum solution using the Davidon- 
Fletcher-Powell Method with G-olden Section Search Technique. 

Cost reduction and improvement m design over a 
recently executed field design has been illustrated. 
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CHAPTER I 


AUTOMATED OPTIMUM DESIGN OP REINFORCED CONCRETE SILO 

1.1 INTRODUCTION: 

Silo is a common storage structure for storing 
granular materials like food grains. The substantial 
losses suffered in the food output because of improper 
or no storage facilities have warranted the attention 
of our country to construct a large number of medium 
sized economical silos. It is this consciousness at 
the national front which has motivated the present 
work. 

Silos are preferred over other conventional 
storages because 1) they are moisture-proof, vermin-— 
proof and insect-proof, 2) they provide the maximum 
space utilization per tonne of the stored material thereby 
reducing the storage cost, 3) they preserve the quality 
of the sotred material? and 4) they entail no loss of 
the stored material while filling aid emptying 
(m handling). 
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Reinforced concrete is preferred as a construct- 
ional material for silo because 1 ) it is economical, 

2) it is fire-proof, and 3) the quality of food gram 
is better preserved. 

Reinforced concrete silos are constructed m 
(*)* 

different ways v J : 

(l) Construction m hollow blocks, 

(n) Pillar ana panel construction, 

(m) Staves silo construction, 

(iv) Pre-cast ring silos, and 

(v) Cast-m-situ silos. 

The first three methods are suitable for relatively 
small capacity. Pre-cast ring silos are used for all 
sizes. Cast-m-situ silos are particularly suitable for 
medium and large sizes and lave been considered m the 
present work. 

The design for a R.C. Silo is essentially a 
trial and error process where the shape, dimension and 
reinforcement details are assumed and the stresses are 
checked for a safe design. Obviously this process 
cannot be repeated many times by hand computation 
because of the sheer labour involved, thus resulting 


* lumbers within such brackets indicate the references. 
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many times m an uneconomical over-designed s true -cure. 

Any saving achieved on the cost of a R.C. Silo is 
going to build-up because of the repetitive requirements 
of this class of structure. Hence, the problem is most 
suitable as an automated optimum design problem. Further- 
more, the forces acting on 9 silo e.g., the dead load, 
the wind load, the earthquake force, and the gram 
pressures depend directly on the shape and dimensions 
of the structure. Thus any attempt of optimal design 
towards the saving of the material cost couples itself 
with the reduction of forces acting on the structure. 

1.2 SCOPE OP THE PRESENT WORK: 

Here, a circular R.C. Silo with conical (self 
cleaning) bottom supported on a system of ring-girder and 
columns as shown m Pige.1 and 2 ha 3 been cast as a 
mathematical programming problem with minimum cost of 
materials used as the objective function. The adopted 
design procedure is the one recommended by Indian 
Standards Specifications^ nJ1< ^ 

The loads on a silo hnve been summarized in 
the laght of IS recommendations^ an< ^ with some 
emphasis on pressures m silos as suggested by different 
investigators, m Chapter II. 
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In Chapter III the analysis and design 
procedure for R.C. Silo has been presented. 

The optimal design formulation for the problem, 
which turns out to have 15 design variables and 61 
constramsts with cost of materials used for the silo 
super structure as objective (or merit) function, is 
described m Chapter IV. 

Chapter V discusses in short the optimization 
seeking methods used for the problem. The Sequential 
Unconstrained Minimization Technique developed by 
Piacco and McCormic giving all the points m the feasible 
region and thus being of great engineering a dvantage has 
been used for converting a constrained problem into an 
unconstrained one. The Davidon-Pletcher-Powell method 
with Golden Section Search Technique has been used to 
seek the solution of the unconstrained minimization 
problem. 

The numerical results have been presented m 
Chapter VI. Reduction m cost and improvement m design 
over an existing field design has been shown. Conclusions 
and scope for future work has also been presented m 
this chapter. 
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In any effort of this hind, the development of 
a computer programme takes the bulk of the effort. The 
automated optimum design programme of a silo developed 
m the present work m fortran language is attached as 
appendix. 



CHAPTER II 


LOADS Off TH3 SILO 


2.1 INTRODUCTION: 

Silo is a deep bin, Bin is a generic term used 
for the class of storage structures e.g., bunker or 
shallow bin and silo or deep bin. The classification 
of a bin into bunker or silo is based on the pressure 
exerted by the stored material 0] ^ to its walls. The 
plane of rupture of the stored material (defined latter 
m this chapter) helps this classification. Most 
significant part of the loads acting on a silo is the 

pressure exerted by the stored material. Various invest- 

. , (4# 5, 6, 7, 8, 9 and 10) . .. , (2,11,12 

igators 9 9 9 9 9 and authorities 

and 13 ) 

have looked into this aspect of silo loading which 
makes it a subject m itself with substantial background 
and interesting history. This has been discussed in 
short m Section 2.3 of this chapter. The other forces 
acting on a silo are 'the dead load, wind load, earthquake 
forces and temperature effects (provided the material 
stored undergoes a temperature change). All these are 
described m the subsequent sections m more or less 
details depending on its earlier documentation. 
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2.2 DEAD LOAD Oil SILO: 

The dead loads on silo comprises of the weight 
of 1) the roof; 2) walls (including partitions, if any); 

3) hoppers or floors; and 4) other structural components 
along with the weight of all other permanent constructions 
and fixtures as well as the mechanical equipments fitted 
on the silo for handling of the material. 

2.3 LIVE LOADS ON SILO: 

Live loads on silo include all loads other than 
the dead loads and may consist of the following: 

i) (a) Live loads on the roof; 

(b) Snow load; 

n ) Wind load; 

in) Earthquake load; 

iv) Pressures due to stored material!* &nd 

v) Loads because of the impact and /ibrations 
due to moving equipments (if any) fitted on 
the silo . 

(2 3 ) 

Indian Standards Specif ications v ’ 7 provide 
the details of the above live loads on silo and these 
have been adopted m the present work. 
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2.3.1 LIVE LOADS OK" THE ROOF: 

(a) Live loadA"^ on flat or curved roof with. 

slope upto 10 degrees to which access is provided is taken 
2 

as 150 kg/m with a maximum of total of 375 kgs. live load 
on roof without proper access is taken as 75 kg /cm . If 
slope (6 r ) is greater than 10 degrees (no access is 
provided), then the live load is computed from the 
expression: 

w r = 75 - 2 (e r - 10) (2.1) 

with a minimum of 40 kg/m^* 

2.3.1(b) Snow Load: If the silo is to be constructed 

m a region of snow fall, it is designed for the actual 

load due to snow or for the live load as specified above 

which ever is more. Actual load due to snow will depend 

upon the shape of the roof and its capacity to retain 

snow. The snow load is assumed to be 2.5 kg/m per cm 

(3) 

depth of the snow . In the present work snow loads 
have not been taken into consideration. 

2.3.2 WIND LOAD: 

The effect of wind is calculated on the basis 

(3) 

of basic pressures given m the I.S. Code for entire 
height of the silo and any projection thereof, having 
due regard to the variation m pressure with height and 
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shape factor. The shape factor for circular structures 

13 ) 

is 0.7 .Basic wind pressure can be had from Fig.lA, 

l>age VI. 1.18 of reference 3 for any height. The value, 

of wind pressure on any structure upto the height of 

2 

30 m is 150 kg/Tom for Northern India (except for Jammu 
and Kashmir). 


2.3.3 EARTHQUAKE LOAD : 

Earthquake shocks cause a movement of ground 
on which the structure is situated resulting into the 
vibration of structures located on it. The intensity 
of ground shock at any location depends upon the 
magnitude of earthquake, the distance of the epicentre, 
duration of ground morion and nature of foundation. 

The response of structure subject to such ground shock 
depends on material, form, size and mode of construction 
of the structure. 

Design acceleration has been specified by X*S. 
Code w ' for structures standing on soils which will not 
settle considerably or slide appreciably due to vibrat- 
ions lasting for a few seconds. Seismic coefficient ( a 
and factor (F) (Page VI. 1.1 3, 16 and 17 of reference 3 ) 
which depend upon seismic zones is multiplied by 
another factor p (Table 5, Page VI. 1.1 4 of reference 3) 
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depending on the soil foundation system. The seismic 

zones for India has been shown m Pig. 13 , Page VI. 1.22 

of reference 3 . Por zone III the values of and P 

are 0.04 and 0.3 respectively. Por raft foundation the 

value of factor (3 is 1.0. As food gram storage structures 

are of post-earthquake importance, seismic coefficients 

and factor P as obtained above are further multiplied 

(x) 

by a factor of 1.5 w . Por concrete structures the per- 

(3) 

centage of damping is taken as 5-10 v . 

In the present work the earthquake force has 
been considered as an equivalent static force acting at 
mid-height of the silo wall. The system has been 
idealized as a vertical cantilever with lumped mass of 
Wfp/g at the free end. The equivalent stiffness is given 
by: 

= IISC (3 E 00l I col /H c b (2.2) 

and the natural time period therefore is 

i n = 2 % /V55i (2 ‘ 3) 

Por the calculated value of T n and assumed 
percentage of damping (5 for the present work); the average 
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acceleration (Accl.) is obtained from Fig- 15, Page VI. 1.16 
of reference 3. Finally, the earthquake force turns out 
to be 


Earthquake = U5 ma3S A ° o1 - (2>4) 

2.3.4 PRESSURES DUE TO STORED MATERIALS: 

Early designers thought that bulk materials 
behave like liquid, and designed the storage structures for 
equivalent liquid pressure. Robert (1882 and 1884) 
demonstrated with the ^elp of experiments on models and 
full size silos that the wall pressures donot increase 
linearly with depth, but part of the weight of stored 
material is transferred to the walls by friction and 
consequently, the vertical and horizontal pressures m the 
lower part of the silo are less than what would be m case 
filled with liquid. Furthermore he found that the walls 
are m vertical compression which is certainly not collaborat- 
ed by a liquid pressure on the walls of the container. 

Janssen (1895) and Airy (1897) came with their 
theories based on statical equilibrium of granul&r materials. 
Large amount of experimental work ha s been conducted and 
reported by Prante (1896), Bovey (1906), Jamieson (1904), 
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Lufft (1904) and Pleissner (1906). letchum ( 1909 ) 

summarised the state of art m a hand book^ 1 ^. There 

he mentions that the lateral pressure due to gram on 

bin walls is less than the vertical pressure (0.3 to 0.6 

of the vertical pressure, depending on the gram etc.) and 

the increase m this pressure is very little after a depth 
1 

of 2 t? to 3 times the width or diameter of the bin. He 
further adds that, m case of the moving gram, the press- 
ure is slightly greater than the pressure of gram at rest 
(maximum variation for ordinary condition is only 10 per cent). 
Prom 1909 until 1950 no significant contribution to this 
study is traceable. It was felt that the problem of 
pressure distribution m silos was resolved for once and 
all. However, refinement m the construction materials 
and improvement m structural design techniques with 
reduced factor of safety caused failures of many storage 
structures designed m accordance with the Janssen's or 
Airy’s theory. Theime^^ has exhaustively reported 
failures of R.C. Silos, and has rightly mentioned that 
majority of failures occurred because of operational 
pressures of storage material , which were much higher 
than the static pressures on which the designs were based. 

Such failures oC storage silos all over the world 
were of much concern among the engineers, and this spurred 
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new investigations into the loading (pressure distribution) 
which a granualr material exerts on the storage structures. 

Prante's (1896) early observations that pressures 
varied widely between the initial condition of charging 
(or filling) and the condition of flow and not constant 
as claimed by Janssen, Airy and Ketchum, were confirmed 
experimentally by Takhtamishev (1938-39 and 1941 - USSR), 
Shumsky (1941 - USSR) , Bernstein (1947 - USSR), Reimbert 
(1943 and 1956-Prance), Bergaft (1959- U.K.), Kim (1959-USSR), 
Kovtun and Platonov (1959-USSR), Pieper and Wenzel (1964- 
Uermany), Turitzin (1963-USA), Blanchard and Walker (1966-U.K.) 
etc. 


It is of interest to note that some of the early 
experimental investigations lead to wrong conclusions as 
they were based on unsophisticated instrumentations. More 
refined testing indicates that the flow pressure exceeds 
the initial (or filling) pressure by a factor of 2 to 4 
and m some cases the high peak loads occur at higher points 
also and not only at the base of the silo. 

In past 30 to 40 years because of the extensive 
experimental and theoretical studies, the behaviour of 
granualr materials have been better understood. This lead 
to the urgent need for modifying the Janssan’s (1895) 
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and Airy's (1897) classical or Reimbert brothers (1956) 
moderate theories for silo pressures based on static 
equilibrium and to account for the dynamic effect m 
flow condition. 

Swiss (Installations d ‘ensilage, 1956, no. 160, 

German (Dm 1055, Blatt 6, Berlin, 1964), Russian (CH 302- 

65, Moscow, 1965), Indian (IS :4995-1968) codes and ACI-313 

( 1 3 ) 

report on bin wall design and construction v still use 
the Janssen formula with little modification here and there 
to account for the increase m pressure based on the 
flow concept. 

('ll) (p) 

The German' ' and the Indian v ' codes arbitrarily 
recommend pressure ratios ^ ^ in filling and in emptying 
as 0.5 and 1.0 respectively. They have recommended 
angle of wall friction ((j) r ) equal to K times the angle 
of internal friction ((])) of the stored material. The 
constant value K suggested, is 0.75 and 0.6 m case 
of filling and emptying respectively. This is obviously 
because of the lack of available test data. 

( 12 ) 

The Russian code' ' recommends that, the values 
of <J) , (j) 1 and should be obtained experimentally for 
the materials m use. To account for discharge condition, 
they have recommended seperate factors for over-pressure 



15 


and serviceability. They have also recommended 
seperate multiplying factors for exterior and interior 
bins. Safarian^"^ has recommended similar multiplying 
factors to be used with Janss©nsnd Reimbert pressure 
schemes. He has suggested and calculated only one factor 
m place of the two factors given by Russian code. 

Theimer^^ has recommended for 50 per cent 
and 100 per cent increase for interior and exterior silos 
over Janssen pressure values for some intermediate 
height. for the rest of the upper and lower part he 
recommends linear decrease from suggested overpressure 
to the maximum values suggested by Janssen. 

All the above schemes indirectly account for 

the flow condition by taking some multiplying factor or 

by increasing the pressure ratio. Walker^ gave an 

approximate theory for pressures and arching based on 

experimental cum theoretical work which accounts directly 

for flow conditions. However, Jenike and Johanson 
( 8 ) 

theory^ gives the only correct mathematical analysis 
of the overpressures in flow conditions. They have very 
clearly pointed to the concentrated reaction at the 
transition from active to passive state of stress field 
m silos, to maintain the field equilibrium. 
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Work is m progress both at the experimental 
and on the theoretical front and universally accepted 
design pressure for the silo is yet to be agreed upon. 

For the present work, the pressure scheme as 

( 2 ) 

per I.S.recommendation v has been used as the basis of 
design, so that the results may be compared with the 
available field designs. However, the programme developed 
can easily incorporate any other pressure scheme without 
difficulty. 

2.4- PRESSURE SEHEME USED IN THE PRESENT WORK: 

The IS:4995 code for Resign of R.C. Bins (Silos) 

( 2 ) 

for Bulk Food Gram Storage v ' gives the scheme for 
pressure on the silo wall and hopper bottom under the 
assumptions : 

a) Material stored m the bin does not undergo 
volume change during storage. 

b) The cohesion for the stored material is too 
small as compared with the internal friction. 

These have been considered m the present work. 

2.4.1 PARAMETERS FOR PRESSURE CALCULATION: 

The parameters required for calculation of 
pressure on silo walls are: 1) unit weight (w) , 2) angle 



17 


of internal friction (tf>), 5) angle of friction between 

wall and stored material (f)), 4) pressure ratio ), 
and 5) height to diameter ratio (h/D). 

2.4.1(a) Unit Weight (w) : - The unit weight of the 
stored material depends upon many variable factors such 
as moisture content, particle size and temperature. Conse- 
quently, universally applicable values cannot be given. 

In order to be exact, the actual value should be obtained 
by test results. However, specifications are given for 
various materials m Table 2 of reference 2 which has been 
used m the present work. 

2.4.1(b) Angle of internal friction (<j)): - The angle of 
internal friction also varies and should be obtained as 
in case of the unit weight. However, here too, the same 
table is referred for the present work. 

2.4.1 (c) Wall friction ((j)') : - It is advisable to 
perform experiments for getting actual values of wall 
friction. In the present work, the angle of wall friction 

has been assumed to be 0.75f) end 0.6(f) during filling 

(2) 

and emptying respectively as suggested by I.S. Code . 

2.4.1(d) Pressure ratio (A ): - The ratio of the 
horizontal to vertical pressures used m the present work 
are 0.5 and 1.0 during filling and emptying respectively. 
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(p) 

The I.S. Code v has purposely suggested these values 
on higher side to account for the over-pressure which 
exists m the real situation (specially m flow condition). 

2.4.1(e) Height to Diameter Ratio: - A silo (deep bin) 
is a tall structure m which the plane of rupture meets 
the opposite wall before it emerges out from the top of the 
filling. The plane of rupture is that surface down which 
a wedge of material bounded by one wall face, the free 
surface and the plane of rupture would start to slide, 
if the boundary walls were to move. The plane of rupture 
distinguishes the shallow and deep bins (Rigs. 3a and 3b) 
and is vital for pressure calculation on the walls and bottom 
of bins. 

In silos, the entire weight of the material 
stored is not transferred to the bottom of the structure. 

A considerable part of this load is resisted by the friction 
between the material stored and the walls of the storage 
structure. This results m the reduction of lateral 
pressure. After certain depth (2^ to 3 times the width 
or diameter of the silo), the increase m pressure with 
depth is very nominal. A silo is preferred to a shallow 
bin because of the above fact. 
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Different authorities have suggested different 

limits on the H/D ratio for deep bin consideration. The 

_ ( 2 ) 

I*S. recommendation v for classifying a bin as a silo 
is as follows: 


H > D tan 0 


(2.5) 


where tan 0 = 



In case H/D ratio does not conform to the above 
equation, the pressure scheme for design consideration 
shall be taken for a shallow bin. In the present work the 
latter eventuality does not arise. 


2.4.2 PRESSURE CALCUL ITI0M5 : 

The maximum values of the horizontal pressure 
on the wall (P^), the vertical pressure on the horizontal 
cross-section of stored material (P v ) and the vertical load 
transferred to the wall per unit area due to friction 
(P ) shown m Fig. 4(a) are calculated using the expressions 
given below m Table 2.1. 
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TABLE 2.1 


(2) 


MAXIMUM PRESSURES 


Name of Pressure 

Luring Pilling 

During Emptying 

Maximum P v 

w.R /(A f ) 

w.R/ (A e >* e ) 

Maximum P^ 


w.r/A* e 

Maximum P 

w.R 

w.R 

w 




2.4.2(a) Variation of Pressure With Depth m Silo: The 
variation of pressures P v , P^ and P w during filling and 
emptying is obtained from the expression v . 

o- e - z/ %) (2-V) 

where, P stands for pressure and suffix l stands for 

v, h or w corresponding to the pressure to be calculated 

and 

= R/(A f A*' ^ ) during filling jj 
= R/(A e /^ e ) during emptying 5 


^i \nax 


is taken from Table 2.1. 



21 


2.4.2(b) Vertical Load on Walls Due to Friction:- The 
vertical load carried by the wall due to friction between 
material stored and bhe wall at any level is equal to the 
product of the vertical load on the wall (P ) times the 
area on which it acts, above that level. It is same as 
the product of total lateral pressure on bhe wall above 
that level (under consideration), and the coefficient 
of wall friction ( ) . 

2.4.2(c) Lateral Pressure hear the Bottom of a Silo: To 
take advantage of the decrease in lateral pressure caused 
by He presence 0 f the bottom of a silo, tho horizontal 
pressure due to emptying, over a height of 1 . 2D or 0.75h, 
which ever is smaller, is taken to be varying linearly 
from the emptying pressure at this height to the filling 
pressure at bottom as shown m Fig. 4b ( i). 

2.4.2(d) Pressure on Hopper Bottom:- The intensity of 

normal pressure ( P ) at any point on the hopper bottom is 

( 2 ) 

taken v ' as 

P^ = P Tr Cos 2 a + p. Sm 2 a + W 0 Cos a (2.9) 

XX v XI o 

where and P^ are obtained from eqn. (2.7) with the 

values of R corresponding to the horizontal section under 
consideration (Fig. 4.b n). 
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2.5 LOAD COMBINATIONS DOR DESIG-IT CONSIDERATIONS : 

To ensure safety and economy m the design of 
any structure, a judicious combination of the working 
loads is necessary, keeping m view, the probability 
of their acting together, their disposition m relation 
to other loads and the severity of stresses or deformations. 

In the present work, the combination of loads on silo 
for design purposes has been taken as follows 

1) Dead load alone 

ii) Dead load + Live load due to storage 
111) Dead load + Lateral load 
iv) Daad load + Lateral Load + live load due to storage. 

The lateral load means wind or seismic force which 
ever is severe and has been considered as the equivalent 
static load acting at the mid-height of the silo. 



CHAPTER III 


ANALYSIS AND DESIGN OP SILO SUPER STRUCTURE 

3.1 INTRODUCTION: 

The silo superstructure consists of the following 
elements : 

1 ) Roof, 

ii) Vertical wall or silo wall, 
ill) Hopper bottom, 
iv ) Ring-girder, and 
v) Supporting columns. 

The procedure for the analysis and design of each 
of the above element is described in more or less detail 
depending on earlier documentation and for the sake of 
completeness of presentation. 

The material for construction considered m the 

present work is reinforced concrete consisting of a 

concrete mix (M-200) and reinforcing mild steel. It may 

( 2 ) 

not be out of place to mention here that I.S. Code 
specifies to use a concrete mix not leaner than M-200. 

3.2 DIMENSIONS OP ELEMENTS : 

3.2.1 DIMENSIONS POR SILO ROOF 

The roof of the silo is designed as a flat 
conical dome. The roof thickness is determined to safely 
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carry the total load (including self-weight of slab, 
waterproofing material and live load) on it. The meridional 
force T^ per unit length is given as 

T 1 = w r h r /(2 sm 2 0) (3.1) 

and the hoop thrust P er length is given as 

T 2 = w r h r Cot 2 0 (3.2) 

Roof thickness (T ) calculated on the above basis, m 
general, comes to be too small and thus; is governed 
by minimum requirement dictated by practical considerations. 
In the present work, a fixed thickness of 8 cm has been 
considered. An additional thickness of 2 cm has been 
accounted for water-proofing material. The slope of the 
roof is taken as 5° and the diameter of the central opening 
as 30 cm. 

3.2.2 DIMENSIONS EOR SILO-WALL: 

The overall size of a silo is governed by H/D 
ratio which depends upon me storage capacity and the 
material stored. As the increase m pressure with depth 
is negligible m the lower part of silo wall, the 
I.S. Coie^ 2 ^ recommends to keep this ratio preferably 
more than two. In tthe present work, the parameter (H/D) 



25 


is a subject of study. The lower limit on it has been 
taken as 1.5 to account for the storage of wheat for 
which the angle of internal friction ((j>) is 28°and 
corresponding minimum H/D ratio cones to be 1.4575 and 
1.4995 m case of filling or emptying respectively. 

For a given storage , tins ratio is automatically fixed 
which corresponds to the minimum cost. The internal 
diameter of the silo (F), has been considered as a design 
variable m me present work, the height (H) gets fixed-up 
for a given storage corresponding to the minimum cost 
design and is calculated as follovrs (Fig. 2):- 

H = H = (total volume of storage - volume stored in 
hopper)/ (it F 2 /4) + H clear (3.5) 

where 

H clear = (3D tan 

The thickness of silo wall (T ) is designed to 

safely carry all the loads coming on it, the procedure 

for which is discussed m latter sections. However, as 

(o') 

specified by the 1-S* Code v , I should not be less than 
the minimum value, l min ( CD1 )> obtained from 

T „ = 15 for D < 6 m (3.5) 

min. — 


= 15 + (D-6)/l .2 for D > 6 m 


(3.6) 
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3.2.3 DIMENSIONS FOR HOPPER BOTTOM: 

The hopper bottom for circular silos is m the 
shape of frustum of a cone . The diameter of the outlet 
opening (D O p 0n ) depends upon the desired rate of withdrawl 
of stored material. In the present worl-c, a constant 
outlet diameter of 30 cm has been considered. In order 
that the silo be self cleaning, the hopper angle, a, with 
the horizontal is to be atleast 10° more'' ' than the angle 
repose of the stored material. For dry and loose fill, 
the angle oC repose and the angle of internal friction 
are approximately equal. In the present work, angle a 
has been taken as 45°. This is because the angle of 
internal friction for materials like wheat to paddy vary 
from 28° to 36°. Once the outlet diameter (D en ) and 
the internal diameter of silo (D ) are known, the 
height of the hopper bottom (H^p) is obtained from 

H hop = tan “-(^cpef/ 2 - 0 (3 ' 7) 

The hopper thickness is designed to safely carry 
the loads coming on it. The design procedure is 
discussed m latter sections. The minimum limit on it 


is taken as 15 cm 
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3.2.4 DIMENSIONS FOR RING- GIRDER : 

The height (H ) and width (B r ) of the ring-girder 
supporting the silo and resting generally on equally 
spaced columns (six in the present work) is designed 
to safely carry all the loads coming on it from silo as 
well as for the interacting forces between ring-girder 
and columns. In the present work, a lower limit of 30 cm 
on these dimensions has been imposed to accommodate the 
reinforcements. Because of practical consideration and 
to use a design table for torsional properties of the 
section, the ratio of these two dimensions has been taken 
m between 1 and 4. 

The analysis and design procedure adopted for the 
ring-girder as presented m the latter sections is meant 
for the idealized equivalent section (Fig. 2). 

3.2.5 DIMENSIONS DOR SUPPORTING COLUMNS: 

The column height (H ) is fixed on the basis of 
clearance required belov* the hopper outlet and is 
obtained as follows (Pig. 2): 


H c = 


V 


cl 


+ B 


hop 


- H 


(3.8) 


The column diameter is designed to safely carry 
all the loads coming on it through the ring-girder. 
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The lower limit on the column diameter is taken as 4-0 cm 
This value keeps the ratio of effective length to the 
lateral dimension of column within 15 (the short column 
requirement, Section 3,8) for an effective length of the 
column upto 6 m. 

3.3 GENERAL DESIGN CONSIDERATIONS: 

The elastic design procedure of reinforced concrete 
has been used m the present work. Stresses due to drying 
shrinkage have been considered for members subject to 
direct tension, viz. , vertical wall and hopper bottom. All 
the elements have been designed on no crack basis (except 
for the ring-girder with large eccentricity). Sufficient 
reinforcement has been provided for members m direct 
tension, designing the same with crack considerations to 
ensure safety even m the worst case. The stability of 
vertical wall has been checked after designing it on the 
basis of strength. 

The values of permissible stresses and modulus of 
elasticity for the constructional materials used m the 
present work have been also summarised m this section. 

3.3.1 STRESSES DUE TO SHRINKAGE * 

( 2 ) 

The shrinkage coefficient , 8 . is taken as 0.0003 

s 

to account for shrinkage effect and accordingly the permissible 
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( 2 ) 

tensile stress m concrete is increased by 33.33 per cent 
when shrinkage stresses are accounted for. 

3.3.2 DESIGN OP A SECTION UNDER DIRECT TENSION: 

The silo wall is subject to hoop tension due to 

lateral pressure (P^). The hopper bottom is subject 
o 

to meridmal and hoop tension caused by pressures P^ and P v » 
the self weight of the hopper and the gram stored m the 
hopper. These elements are designed on no crack basis. 
However, sufficient reinforcement is provided at each 
cross-section so that, if a section cracks due to unexpected 
operational pressure or lack of supervision at the time of 
construction, the total tension can be taken by the 
reinforcing steel only. 

If T is the tension per metre length, then the 
thickness of cross-section (t) is computed as mentioned 
below: 

Actual tensile stress m concrete ( £, ± ) is given as 

Ci = T/(A C + m A g ) <^ t (3.9) 

Substituting, T = f . A_ and A = 100 t 

f s A s < (100 t + m A g ) 
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or . 


\im ^s f s ^s ” t )/( 100 \) 

= T (f s - m <^)/(100 f s <^ t ) (3.10) 


The expression (3.10) is known as Portland cement 
formula C 1 5 » 1 6 ) 

xuxmuj-d and ls derived on no crack basis. The 

thickness provided should be atleast equal to t or the 

mm 

minimum thickness specified m section 3.2. 

When shrinkage is also accounted the equations 
(3.9 and 3.10) get modified to 


= (T + 6 S E g A )/ (A + m A _) 


3.11) 


and \m = T -s 


(f e + e s E s " m ^1)/(100 + ) (3.12) 


S v t‘ 


3.3.3 STABILITY CONSIDERATION FOR SILO WALL: 

In general, silo is a slender structure. In 

addition to strength considerations, check for stability 

for vertical wall under different loading conditions is 

necessary to ensure the safety against buckling failure. 

Stability of silo wall has been checked m the present 

work by limiting the compressive stresses as specified 
( 2 ) 

by I.S. Code v J . The compressive stresses are calculated 
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on the basis of net cross-section (deducting all 

openings and properly accounting for reinforcement). The 

total compressive stress calculated on the above basis 

with lateral load and without lateral load should not 

exceed 0.2 f and 0.15 f. respectively, 
c c 

3.3.4- PERMISSIBLE STRESSES A ED MODULUS OE ELASTICITY: 

The permissible stresses for reinforcing mild 

steel and concrete (M-200) as used m the present work 

(2 S') 

have been taken from the I.S. Code v * ' and are as follows: 


Particulars 


Rotation 


Permissible 

stress 


<p 

(kg/cm - ) 


l) General Reinforcements : 


Ordinary reinforcement 

f 

o 

1300 

(other than helical) 



Helical reinforcement 

f sh 

1000 

Shear reinforcement 

f ss 

1250 


ii ) Silo Reinforcements 

Reinforcement m direct tension 0 

o 

Reinforcement m bending tension x 
on exposed face 5 

Reinforcement for members less x 
than 22.5 cm m thickness 2 


1000 
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Particulars 

notation 

Permissible 

stress 



(kg/cm 2 ) 

ii) Reinforcement for meml ers 22.5 cm 
or more m thickness (other than 

f s2 

1250 

those specified above) 



lii) Concrete (M200) 



Cube strength of concrete at 

28 days 

f c 

200 

Bending compression 

<f c b 

70 

Direct compression 

<Tc 

50 

Bending tension 

tftb 

17 

Direct tension 


12 

Shear 

& 

17 


The modulus of elasticity of steel (E_) is taken as 

s 

C' o 

2.1 x 10° kg/crn . The modular ratio (m) is obtained as 
m = 2800/f (3.13) 

O 

and the modulus of elasticity of concrete (E c ) Is obtained as 

E„ = E /m 
c s 


(3.14) 
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3.4 ANALYSIS AND DESIGN OE SILO WALL: 

The wall thickness and hoop reinforcement are 
designed for hoop force (E^ es ) developed due to the lateral 
pressure (P^) and ^s given as 

P de s = ^h 11 / 2 (3.15) 

where , P h is the absolute maximum pressure m either 
case of filling or emptying and for any point on the 
wall . 

The wall thickness (T w ) is to be designed m 
such a way that it satisfies the no crack requirement 
as well as the minimum code requirement as discussed 
earlier m sections (3.3.2 and 3.2.2). 

The hoop re inforce me nt (A g ^ ) per meter length is 
calculated treating that the hoop force is taken by steel 
only and is obtained as follows : 

■^sl " ^des^sl (3.16) 

The hoop reinforcement as calculated above 

should be more than the minimum requirement, i.e., 

( 2 ) 

0.3 per cent v ' of the cross-section. 

The verti-cal reinforcement (A^) in "the silo wall 

(p) 

with a minimum of 0.3 per cent v J is provided, such that 
the section is sufficient to resist the forces coming on 
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it m either case, namely l) empty silo with lateral load, 
2) silo full with no lateral load, and 3) silo full with 
lateral load. Since the bottom most section of the silo 
wall experiences the maximum stresses, it is this section 
which is checked. No curtailment of vertical reinforcement 
is proposed, fimce^the reinforcement at this root section 
generally corresponds to near minimum requirement. 

The stresses at the root cross-section are 
calculated as follows: 

The compressive stress m concrete due to self 
weight of the wall and the load on it coming from the 
roof is given by 

^sw ~ < ¥ top + V, T w w c^ // ^‘equ (3.17) 


where, A represents the equivalent area per meter 
^ 

length of the perimeter of the ring and is given as: 


A 


equ 


= T + (m-1 ) L 
w a 


s2 


(3.18) 


The compressive stress m concrete due to the frictional 
load coming on the wall due to stored material is given by : 


S sg ~ T ff/ A equ (3.19) 

where, T^ is the total maximum frictional load on 
unit penphenal length, of silo wall and is given by: 


T ff = < h tw - H clear> * ' P v 


( 3 . 20 ) 
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where, P T is the vertical pressure on the horizontal 
cross-section of the stored material at the junction 
point of the vertical wall and hopper. 

The lateral load (hj 10rz ), maximum of the equivalent 
static earthquake force and the wind force, is considered 
to be acting at mid height of the silo wall. Therefore, 
the maximum bending stress is obtained as: 


S sh = ^horz. ’hP/'W 


( 3 . 21 ) 


where, I u represents the approximate moment of inertia 
of the equivalent ring and is given by 


3 " 

equ 


D 3 (T + (m-1 ) ) /8 


The combination of these stresses are considered 
to correspond to the abo^e said loading conditions and 
accordingly the following cases are checked to ensure that 
the stresses are wibhm the limit : 


I) Empty silo vvith lateral load: For no crack the 
relation 

S sh - S sw < ifx < 3 ' 22) 

is to be satisfied. 

II ) Silo full with no lateral load: The stability 
consideration resulting into the reduced stress check as 
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per section 3.3 is satisfied by: 

S sw +S sg < 0.15 f c (3.23) 

ill) Silo full with lateral load: In this case also 

the stability consideration as per section 3.3 is satisfied 

by: 


sw 


S sg + S sh 


< 0.2 f 


(3.24) 


3.5 ANALYSIS AND LESION OE HOPPER BOTTOM: 

The conical hopper bottom is designed to safely 
with- stand the meridional force (T ^ ) and hoop force (Tg). 

Total meridional tension (T^) per unit length of 
the periphery of the hopper bottom is obtained by consider- 
ing the vertical equilibrium of forces at that level 
(Fig. 5a) and is given as: 

1,51.. = .l| V* + V + \ < 5 - 25) 

In the present work the meridional force (EM^ ) at the 

o 

junction of the hopper and vertical wall and the meridional 

force (EM 2 ) at mid hopp#x height are calculated. The 

maximum of these twr t os per unit length is taken as the 

0 

design force (3H% es )' The mendinal reinforcement and the 
thickness of hopper bottom are computed exactly the same way 
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as the case of vertical wall subject to P^ gg . 

Referring to Rig. 5(b), the hoop force (Tg) 
per unit length is obtained from the equilibrium equation: 

T 1 /R 1 + T 2 /R 2 = P n (3.26) 

For a flat conical hopper, R^ is infinite and R 2 is 
I>h/(2 Sm a). So the hoop force is given by 

T 2 = P n R 2 = (P n \ )/(2 Sin a) (3.27) 

In the present work, the hoop force at the junction 
of the hopper and vertical wall and at mid hopper height are 
calculated. The maximum of these two forces per unit length 
is taken as the design force (FH^gg)* The reinforcement is 
computed and the thickness is checked m the manner 
similar to one desoribed earlier. 

The meridinal and hoop reinforcements provided 

m the hopper bottom should not be less than the minimum 

( 2 ) 

requirements, i.e., 0.3 per cent of the cross-section^ . 

3.6 EVALUATION OP INTER- ACTING FORCES RETWEEN RING GIRDER 

AND COLUIfNS : 

The analysis of the sub-structure consisting of the 
circular ring-girder and equally spaced supporting columns, 
forming a staging or supporting structure for the silo 
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( 1 7 ) 

has been presented on the basis of Safarxan’s work'- . 
Assumptions have been made to simplify the analysis, 
viz. , 1 ) the pentagonal section of the ring-girder is 
replaced by an equivalent rectangular cross-section 
having the same area and the same centroid (Pig. 6); and 
2) the centroidal axes of the columns are assumed to pass 
through the centroid of the equivalent section of the 
rmg-grider and with the mean diameter circumference of 
the silo wall cross-section. In otherwords the silo 
wall, the ring-girder and the columns are not placed 
eccentric to each other. 

3.6,1 EQUIP ALE NT SECTION OF THE RING— GIRDER : 

The gross cross-sectional area 'of the ring-girder 
(A r ) is given by 

A r = a-j b^ - a 2 b 2 /2.0 (3.28) 

The coordinates of tne centroid of the pentagon 
from the origin 0 (Pig. 6) are given as 

X 1 = (a 1 b^/2 - ( a? b 2 /2) (b 1 - b 2 /3)J /A r 

(3.29) 

and = |a^ b^/2 - (a 2 b 2 /2) (a^a^)] A r (3.30) 
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The height (H ) and width (B r ) of the equivalent 
rectangular cross-section is given by 



ffi 

i! 

$ 

-A 

(3.31) 

and 

B r = V K r 

(3.32) 


3.6.2 LOADS OF THE RING-GIRDER AND C ODUMS: 

As shown m Fig. 6, the force S which is 

o 

equal to the meridional tension, EM^ , coming from 
hopper to the ring-girder, can be resolved m horizontal 
and vertical components as S' and S'. The net unit 
horizontal force (S x ), trc total unit vertical load (Sy) 
and the uniform torsional moment (M^) shown m Fig. 6 
and acting on the ring-girder are obtained as follows: 



S x 

- s 'x 

- p h 

(3.33) 


S Y 

= S'y 

+ < W top + \all + T ff + V 

(3.34) 



= w T / 

( 2tc r m } 

(3.35) 

and 

% 

= S e 


(3.36) 

where , 

e 

II 

H 

+ 

T w Sin a - T h )/2 

(3.37) 


The basic supporting frame with loads coming 
on it is shown m Fig. 7(a). The free-body sketches 
showing inter- acting forces on ring-girder and column 
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are given m Pig. 7(b). These inter-aeting forces are 
obtained by satisfying compatibility as follows: 


1 ) Deflection of rmg-grider = Deflection of column 


or, 


A 


A 


+ 


A 


A 



(3.38) 


y 

where ^A 



and 



is ihe deflection of ring due to S„(always positive), 

.A. 

is the deflection of ring due to (positive 
for the direction shown m Pig. 7b). 

is the deflection of column due to H. and 

A 

M (positive inward). 


li ) Rotation of rmg-girder = Rotation of column 


or 9 


A 


+ 9 


= 0 


* » i 


A 


(3.39) 


where , 0' A is the rotation of rmg-girder due to the 
applied moment M+ (always positive), 


9 » ' 


is the rotation of rmg-girder due to 
(negative for M. as shown m Pig. 7 }, 


and 9 ,,f » is the rotation of column due to H ; and IL, 
A A -a 

(positive for as shown m Pig. 7b). 

The expressions for the various displacements 

(17) 

and rotations described above are given below : 
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x 


1 ) The uniform radial displacement A . of the 

*** li 

circular ring— girder due to uniform horizontal force S, 
is obtained by hook’s law as: 

(3.40) 


A = 

A 


S x E P 


2 ) 


The radial displacement 



of circular ring-girder 


at the point of force application for any number of equally 
spaced equal forces, 11 is given as: 





K’) / (2 £ r 



(3.41) 


The value of K' is calculated corresponding to the 
value of ^ (Fig. 8) from the expression: 

K’ = + (Sm lj) 1 Cos fy )/2 ) l/Sin 2 ^ - l/fy 

(3.42) 

= 0.003364 for six columns as (jh - a/6. 

The value of K’ for any other set of columns, e.g. , 

4, 8, 10, or 12 are given m Table 3.1. 


3) The angh of rotation 9'^ "'■ p the circular ring- 

girder of uniform rectangular cross-section subject to 
uniformly distribute! twisting moment, M^» along its 
centre line is given as: 
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6 'a IvI t r m^ ^ E r H r 1 n ^ r 2^ r 1 ^ (3.4-3) 

4) The angle of rotation 9 " of the support points of the 

lx. 

ring— girder due to interacting moment , at each support 
is determined using superposition of the -values for 
pairs of equal and diametrically opposite applied moments. 
Fig. 9 shows a ring-girder with one pair of such moments. 

The bending and torsion at any section of this 
ring (Fig. 9) are given as: 

(M b ) p = (M ta Cos <J> 2 )/2.0 (3,44) 

and (Mj) ^ = Sin (j) 2 )/2.0 (3.45) 

The angle of rot ’tion at any section of the ring 
- (Fig. 9) is given as: 

q , = — ^ — - a l — f (J) 0 - Sin 2 ^2 : (3.46) 

^ Sc* L 2 "2 J 

where, 0 < (j)2 _< m/2 

At the either support (j> 2 is m/2 and accordingly the 
angle of rotation is given by 

0 „/2 = (M TA tn °0 


(3.47) 
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Consequently, the general expression for the angle of 
rotation at the support for various numbers of 

diametrically opposite moments is given as 

e'/= (M ta r m x K")/(16 C) (3.48) 


where , K ,f is a numerical factor and its value is 
calculated on the principle of superposition depending 
on the number of supports (Pig. 8). The value of the 
coefficient K' 1 for six columns is calculated below: 


| K” = 



Sm 2 




,%/2 


+ 2 



n/6 

Sin 2^ I 



o 


or, 


n 

2 


K« ' = 


( | - 0) *■ 2 ( g - 


13/2 v 

- T — j 


or, k ’ ’ =1.1153 

The value of K ' ’ for any other set of columns, e.g. 
4, 8, 10 or 12 are given m Table 3.1. 

TABLE 3.1^ 17 ^ NUMERICAL C0EPF1CIENTS K’ AND K** 


No. of supports 

K ' 

K » » 

4 

0.012159 

1 .0000 

6 

0.003364 

1.1153 

8 

0.001387 

1.3633 

10 

0.000701 

1.6199 

12 

0.000404 

1 .8972 
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5) The horizontal displace lent ^ and the angle of 
rotation © ' ' '^ at top o_’ the column are calculated 
treating the same as a cantilever. Therefore, 




and 9 


! I I 


A 


H A H c 


3 E c ol Acl 


H A H o 

2 B eol T ool 


+ 




^A I 
2E col 

j 2 

i£L 

I ool 

(3.4-9) 

Ia 

H c 

(3.50) 


E col I col 


/ // , IJJ 

Substituting the values of /\^, /.V ^ and /j, ^ m 
equation (3.38) ard the values of 9 8 ' and O’’'^ 

m equation (3.39) we get: 


S r 2 
x m 


A E 
r r 


+ 


H r^ 
n A m 

2-^x ^r 


K» =- 


H A H c 


+ 


^A 


h; 


3E col I col 


J col col 


(3.51) 


. 12 M, r 

and t m 

E r H r ln ( r 2^ r 1 ^ 


ii m . r tc 
+ Tfr^-r K " 


H. 

A c 

^ E col W 


+ 


> H ° 

E col T C0l 


(3.52) 


Putting the value of torsional rigidity, c‘, as 
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c» = & c J r = ^°* 425 * V ( a 1 B r } 


(3.55) 


and taking same value for modulus of elasticity of concrete 
used for the ring-girder and for columns (i.e. E r = E c0 ^)t 
the above equations simplify to 


( 



K 1 + 





21 


col 



s x 4 



(3.54) 


h: 


and 


H 


21 


A 


+ ( 


r m 71 


col 


6.8 a., B + 


I" - 


K 12 M, r 

) m-, = - - t B — 

I col ^ ln(r 2 /r 1 ) 


(3.55) 


where, a. is the torsional coefficient for the rectangular 
cross-section and is given m Table 3.2. 

The solution of the simultaneous equations (3.54) and 
(3.55) gives ana The area and moment of inertia 

of the reinforced concrete cross-section of the ring-girder 
and the columns used m the above equations are computed 
m the present work on the consideration of the overall 
cross-section ignoring the reinforcement. 
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TABLE 3.2 


TORSIOMh PROPERTIES OB RECTaESUIuIR CROSS-SECTION 


(17) 


Section 


Torsional 
const an! 


Torsional Values of coefficients 

section — — ^ — — — — 

modulus and Ratio = _r a p. 
maximum E 

shear stress r 


S T =p B l 

1 .0 

0.140 

0.208 


1.5 

0.294 

0.346 

M, 

tors 

2.0 

0.457 

0.493 

1 ’ = S m 

max T 

3.0 

0.790 

0.801 


4.0 

1.123 

1.150 


3.7 ANALYSIS AND EES IGF OP RING— G-IRBER : 

The ring-girder is a doubly reinforced curved 
beam. The scheme for the evaluation of thrust, shear force, 
bending moment and twisting moment at particular sections 
of the ring-girder, e.g., support point, mid-span point and pom" 
of maximum torsion has been presented. A design table has 
been prepared for the ease m computation of these forces 
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and moments. This table can be used as a supplement to 
Table 3 of reference 2. The design procedure for a 
girder subject to shear and torsion has been presented. 
Formula has been derived for a doubly reinforced girder 
subject to bending and axial thrust m terms of non-dimen- 
sionalized parameters for the calculation of stresses 
both for small and large eccentricities. The stresses 
developed due to forces acting on the girder have been 
checked against the permissible stresses. The reinforce- 
ments provided m the girder have been checked against 
the sum of the reinforcements required on the basis of 
balanced section to resist the forces and moments acting on 
the girder. 

3.7.1 EVALUATION OF FORCES AND ''0~ir T ~3 : 

i) The total load (W^) or the uniform vertical 
load (Sy) acting on the ring-girder causes bending perpendi- 
cular to the plane of the ring-girder and a torsional moment 

on it. These moments have been determined on the basis of 

(2 17) 

the coefficients 1 which are obtained from the Ketchum's 
equations^ 

li) The uniform horizontal force (S ) causes a 

circumferential thrust which is constant at all cross-section 

and is equal to S v r . 

x hi 
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in) The uniform twisting moment (M^ ) acting 
on the ring-girder due to the eccentricity of the meridinal 
force transmitted from the hopper to the ring-girder causes 
bending moment perpendicular to the plane of ring-girder 
and is equal to M^. r . 

iv) The concentrated moment (1^,^) acting at the 
supports subject the ring-girder to bending perpendicular 
to the plane of the ring-girder and to a torsional 
moment on it. These values at any cross-section are 
obtained using the principle of superposition from 
equations (3.44) and (3.45). 

v) The horizontal force (H^) acting at the supports 
subject the rmg-girder to a bending m the plane of the 
ring and to a non-uniform circumferential compression. 

In general, stresses due to the lorce (H^) are relatively 

(' 17 ') 

small v and have been neglected m the present work. 

The coefficients for vertical shear force, thrust, 
bending and twisting moment at the supports, mid-spans and 
at the points of maximum torsion, for a ring-girder 
supported on 4, 6, 8, 10, or 12 columns have been 

( 17 ) 

calculated on the basis of the values given by Safanan 
and are presented m Table 3.3. 



TABLS 3.3** SCHEMATIC TABLE FOR THRUST; SBBAREORCE . BENDING MOMENT AND TORQ.IJE CALCULATION AT VARIOUS 
; : . AT KEY POINTS IN A SILO RING-GipER ■ ' 


No. of 
suppprt- 
ing 
cols, 

NSC 

Vertical 
load on 
col. 

Vertical £ 

3.E. due to- 11 

Bending Moments 




\ Torsional moment 

Angular 
distance from 
• columns to 
the point 

< P mJ* ° f 

maximum 

Torsion 

at cols. 

k,w t 

•f'at the point 
«>*nt> 

; "2 \ 

_Due to W T 


Due to yt 


Due to W T 

Due t.OiM^ 

at col, 

^^m 

at midway 

between 

col. 

K 4 W T V 

at col, 

■’ ®s “ta 

at midway 

between 

cols. 

^a 

at the point at the 
(P +)* ? oin N 

4"ta SgV 1 ,. 

at tiie 
point . 

*9 %A 



E 1 

^2 


K 4 

*5 

K 6 

4 

*8 ' 

_ V . 


1 

2 

. - 3 

. : v • 4" 

-5 ■ 

6 

7 

. 8 

9 

10 

11 

12 

4 

V 4 

.125 

.07167 

'-*05415 

0,01762 

-.5000 

-*7071 

-.6366 

*0053 

*3069 

19 ° 12 V 

• 6 


,08555 

*04797 

-*01482 

O.OO 751 

-* 8660. 

-1.0000 

-*9549 

.00151 

.2982 

12°44' 

8 

V 8 

» 06250 

*05597 

-*00827 

0.00416 

-1.2071 

- 1*3071 

-1 .2880 

.00064 

.2928 

9°33 ’ 

10 

Wrp/10 

*05000 

*02884 

-*00.527 

0*00266 

-1.5390 

- 1 . 6.170 

-1.5915 

.00032 

*2916 

7° 37' 

12 

Wrn/12 

.04167 

.02405 

-.OO 565 

0.00189 

- 1.8660 

- 1.9313 

-1.9049 

.000185 

.2505 

6°21 * 


Note- On the ring-girder - the compressive force due to S is S r On the ring, girder - the vertical shear force dur to W^at the mid 

x x m point between columns is zero, 

- the bending moment due to- is M t r^ _ the bending moment due .to W y at points (P mt )*is zero 

— the torsional moment due to at the support and mid— span _ torsional moment due to ML, , at columns and. mid 

' points are zero. span points are 0 . 51 ^, and zero iA respectively. 


DESIGN FORGES & MOMENTS; . . ' 

force = S 3^ Vertical shear force /= K, Wf at the support points 
X > = ^ at the points (P^)* 


Compressive 
Torsional Moment 


/ *8 W T r m + *3 “ta 


* Refer to Table column no... 12. 


°.5% 


V 

at the points (p )* 

at the columns (or support points) 


Bending moments 


r 


“3 m 3 "M rt -m- points 

K 4 Vm + Vu + V. at the 

*7 % + “t r m at 

v mt-' 


** coefficients for this Table has been calculated on the basis of Safarian's work {l7 > to match with the I.S. Code< 2 > 
h ' h’ ani % Dre same 53 " tJlG coefficients re . K a t Table 3 01 Ref* 2. 


Here, the coefficient; 
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3.7.2 DESIGH OF GIRDER FOR SHEAR AND TORSION: 

Shear stresses (q ana q * ) m the girder due to 
shear force and torsional moment are respectively 
calculated as 


1 = Q/(P r Dd) 

and 

q' = M t / (p 1 p^) 


(3.56) 

(3.57) 


where, the coefficient p^ is to he taken from Table 3.2. 

If (q + q f ) is greater than the permissible shear stress 
(<)- shear reinforcements (stirrups) are provided, 
otherwise not. However, if (q + q') exceeds four times the 
permissible shear stress m concrete, the section is 
redesigned. 

Reinforcement for shear: If A w1 be the area of 

the reinforcement for shear, the shear resistance developed 

(3) 

is given as : 


Q* = £ss Aw.j qd/p (3.58) 

and consequently the required area of stirrup per unit 
length for shear force (Q) is obtained as: 


Aw 1 

_ 



Q 

f ss 3d 


(3.59) 
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Reinforcement for torsion:- The area of stirrups 
(Aw 2 ) required for a torsional moment (M^) is given^^ by: 

Aw p 

= Aw g = %/(0.8 f sg x y ) (3.60) 

In addition to the shear stirrups (Awp) longitudinal 

f 3 ) 

, reinforcements are also provided^ to take care of the 
diagonal tension. These reinforcements have the same 
volume per unit length as contained m the corresponding legs 
of the stirrups and are provided along with the longitu- 
dinal reinforcements m addition to those required for 
bending. 


The total stirrup reinforcement provided per unit 
length is 

Aw = Aw^ + Aw 2 

~ Ar ^sr3r (3.61) 

3.7.3 DESIGN OF GIRDER FOR BENDING l HD AXIAL THRUST: 

If P ind M be the axial thrust and bending moment 
respectively the eccentricity is given by 

e x = M/P 


(3.62) 
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This eccentricity is temed as small or large depending 
on whether it is smaller or greater than one sixth times 
the depth of the girder. 

In case of small eccentricity complete section 
is supposed to be effective and the design is carried out 
on no crack basis. The centre of gravity line for the 
section is calculated by equating the moment of the 
equivalent area on both the sides of this line (Fig. 10a). 

B r n2 /2 + (m-1) A c (n-Cov) = (H r - n) 2 B r / 2 
+ (m-1) (H -n - Cov) 

Solving this the depth of the neutral axis is obtained 
as : 

Cov (m-1 ) il c + H 2 B r/2 + (m-1 ) ^ (H r - Cov) 

(m-1) A c + H r B r + (m-1) ^ 

(3.63) 


Nou*dimensionalizing both the sides by 
putting N = n^ r> A srl = ^/(.H^) and = V( H r B r>’ 

the above equation (3.63) simplifies to 

0.5 + (in— 1 ) (Ag ri Cov/H r + Asr 2 0-Cov/H r )) 


1.0 + (m-1 ) (h sr -j + Asr 2 ^ 


(3.64) 
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The equivalent area of the girder shown m Pig. 10(a) 
is given by 


A e “ H r B r + ( m ~l) (A c + A^) 


'J 

1 .0 + (m-1 ) (Asr-j + ^r 2 )J H r B r (3* 6 5) 

4 


The moment of inertia of the equivalent section about 
centre of gravity line is given by: 


I e = B r V 1 2 + B r H r ^r/2" n) ' 


+ (m-l) A c (n-Cov)^ + A^ (H r - Coy - n)‘ 


B r H r f Vl2 + (°* 5 " N)2 + (m-1)^A sr1 (F-Cov/H r )‘ 


+ A sr? _ (1~Cov/H r - N)' 


( 3 . 66 ) 


The compressive stress due to axial thrust is same 
throughout and is given as: 


S = b /A 

- c ' e 


(3.67) 


The maximum compressive and tensile stresses 
due 1 to bending me lent are obtained from 


/>cb= M NH r /:C e 


( 3 . 68 ) 
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and 

' r tb 
respectively. 


(1-N) H r /I g 


(3.69) 


The net tensile stress due to combined bending 
and thrust is given by 


S t = M (1-N) K r /I e - P/A e (3.70) 

should not exceed the allowaole tensile stress 
( ^tb ) m order that the concrete does not crack. 
For no compression failure the relation 



6"c 3^ct 


should be satisfied. 


(3.71) 


In case of large eccentricity, excessive tensile 
stress is developed causing cracks m the section and 
the design is carried out with crack considerations. 
Hero, concrete m the tension zone is neglected. The 
Fig. 10(b) shows a girder section subject to a load 
with large eccentricity. The stresses at a point is 
proportional to its distance from the neutral axis, and 
is given ns: 


s 


6 . y 


(3.72) 
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where, y = distance of the point from the neutral axis, 

end 6 = a coefficient of proportionality. 

For equilibrium, 

total compression - total tension = Thrust on the 

section 


or. 


n 

/ B &.y dy + (m-1 ) A_ 5* (n-Cov) 
o 0 


- u Aj 


6 (d - n) =1 (3.73) 


or, 


P = 6 j B r n 2 / 2 + (m-1) A c (n-Cov) -m A^(d-n) 


(3.74) 


Similarly, 

moment of resistance = bending moment 

or, total moment of the induced 5 

compressive and tensile stresses | = bending moment 

about the mitral axis I 

A * 


n 

’ "r 2 


or . _ 6. (2/3 )n + (m-1) L Q d/n-Cov)^ + m A^.(d-n)^6 


= P c 


n 


(3.75) 
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Non-dimensionali zing both, the sides, the equation (3.78) 
simplifies to 


11 (IT + B N 4 G) - A = 0 


or N = A/(h + B N + C) 


(3.79) 


where 


A 


ry 

= 6 /t _1 ^ A sr1 Cov / H r + m A sr2 ( 1 “ c °VH r ) 


(e x /H r - 0.5) 


^ (m-1 ) A sr1 (Cov/H r ) + m A sr2 (1-Cov/^) 



B = 3 (e x /H r - 0.5) = 3 e x /H r - 1.5 


0 = 6 ; (e x /H r - 0.5)) A w1 (m-1) + A ov0 m/ 


t 




srl 


sr2 


+ A sr1 (m-1) Cov/H r + A gr2 m (1-Cov/H r ) J 


The value of H is obtained by iterative technique from 
equation (3.79). 

The total compressive stress and the net tensile 
stress due to combined bending and thrust is computed 
as follows: 
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Talcing moment about the tensile reinforcement we get 

B r n (c/2) (I-I r - Cov-n/3) + (m-1 ) A Q ) c (H^-2 Cov) 

= P ^ e x + H r/2 “ Cov ) 

or e _ p (e x + F -r/2 - c ° ^ _ 

0.5 B r n (H r -C0T-n/3) + (m-1 )k 0 (^=~) (H p -2 Gov) 

P (® X /K r + 0.5 - Cov/H r ) 

B r H r f0.5 N(1-Cov/H r -N/3) + (m-1) A sr1 ' 1-Cov/(lK r )\ (1-2Cov/H 7 

L t ' 

(3.80) 

The value of c gives the actual maximum compressive stress 
on the section and this snoald not exceed the permissible 
compressive stress m concrete 

The maximum tensile stress is calculated from the 
triangle of stress as follows: 

m c _ n _ N 

h; H~-Go \r-n 1_ Cov/H -N 

r ' r 

or t = m c (1-N - Cov/H r ) (3. 81) 

“ U 

This value of t should not exceed the allowable tensile 


o o Q 

O bit DO 


in steel. 
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When the girder is subject to negative or 
hogging moment, the upper re inf orcement becomes tensile 
reinforcement and the lower one becomes compressive 
reinforcement and accordingly the formulae for calculating 
the stresses get modified. 

3.7.4 REINFORCEMENT REQUIREMENT FOR A BALANCED SECTION: 

A balanced or critically reinforced girder 
section is one m which the extreme concrete fibre m 
compression and steel m tension reaches their respective 
permissible stress simultaneously. For a balanced 
section the coeffic^e nts for the depth of the neutral axis 
(N ) and the lever arm ( j-. r ) are given as 

M cr = 1 A 1 + f s /(m ^cb ) (3.82) 

and J cr = 1 - N cr /3 (3.83) 

For a balanced section subject to a moment II, the area 
of tensile steel (A 0 ) required for reinforcing the 

o 

section is obtained as 

A s = ?!/(f s J cr (H r - Cov) ) (3.84) 

As the ring-girder is subject to shear force, 
twistin' 1- moment a nd bending moment, simultaneously, the 
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total reinforcement requirements are checked against 
the reinforcements provided* in all such possible cases. 

The stirrup reinforcement provided at a section should 
be at least equal to the stirrup reinforcement required at 
that section to safely carry the shear force as well as 
the torsion as discussed m Section 3.7.2. The longitudinal 
reinforcement provided at a section should be at least equal 
to the sum of the reinforcement required due to bending as 
well as due to torsion. The required reinforcements are 
calculated with balanced section consideration. 

The lower limit on girder reinforcements, viz., 
upper longitudinal, lower longitudinal and stirrups 
has been taken as 0.3 per cent of the girder cross- 
section. 


3.8 ANALYSIS AND DESIGN OD COITFUIS : 

The silo super structure has been supported on 
a sot of six equally spaced spirally reinforced concrete 
column m the present work. The design procedure for a 
spirally reinforced concrete column is discussed. The 
columns are designed to carry lateral load on the silo with 
and without the stored material m the silo. 
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3.8.1 SILO UNDER LATERAL LOAD : 

As described earlier, the lateral load considered 
m the present work is the maximum of the wind or earth- 
quake force, Ehorz* end acts at the mid-height of the 
silo. The longitudinal force m the column due to 
bending effect caused by the lateral load is considered 
to vary linearly with the distance from the line 
perpendicular to the wind direction or the direction of 
earthquake and passing through the centroid of the 
group of columns. The absolute maximum possible tensile 
or compressive force which may occur m a column is given 

t>y 

F hors. (H r + t V»/2 ) (3.85) 

lodh 

In addition to this all the six columns are subject to 
horizontal force 


E 

h 


= E 


horz /6 


( 3 . 86 ) 


m the direction of wind or the direction of earthquake. 
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3.8.1 (a) Silo Full with lateral load: - In this case 
one of the column is subject to maximum compressive 
force which is given as 


! loai = V 6 - 0 + C lodh 


(3.87) 


The moment acting on the top of the column is given as 


C bm1 “ ^TA 


( 3 . 88 ) 


In addition, every column is subject to horizontal 
forces, and *h m the radial direction and in 
the direction of wind or the direction of earthquake 
respectively as discussed earlier. In a particular case 
these two horizontal forces are added together. Therefore , 
the base of the column is subject to a downward load of 


C loa2 = C lod1 + w c A col H c 
and maximum bemdmg moment of 


(3.89) 


C bm 2 = K H A + V H c - % I 


(3.90) 
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3.8.1(b) Empty Silo with lateral load:- In this case 
one of the column is subject to maximum tensile force 
which is calculated as follows : 

The downward load on the column for empty silo 
is given as : 

^lod3 = “ height of total material 

stored m the silo)/6 

(3.91) 

So, that net tensile force m the column is given by 
C lod t = C lodh “ C lod3 (3.92) 

It may be mentioned here that m the case of empty silo 
the bending moment coming on to the column is m general 
very small and hence is neglected. 

3.8.2 EFFECTIVE LENGTH OF THE COLUMN: 

The effective column length (H c0 ^) is taken 
as ^fact ' t:Lmes ike unsupported length of column (H c ). 
The value of this factor (Elj ac ^-) depends on the end 
conditions of the column and is taken in accordance with 
Table 8, Page VI 5.19 of reference 3. In the present 
work, the end conditions are considered to be fixed at 
the base and partially fixed at the top for which El^ ac ^_ 
has been taken as 0.8. 



64 


3.8.3 SHORT AND LONG COLUMN CONSIDERATIONS: 

If the ratio of effective length of the column 
to the least lateral dimension (radius of gyration) 
is less than 15 (50), the coluiin is termed as short 
column, otherwise, it is termed as long column. 

The permissible stresses m a reinforced 
concrete long column shall not exceed those which 
result from the multiplication of the appropriate 
maximum permissible stresses and the coefficient c r> 
obtained from 



(3.93) 

(3.94) 


where, D mln and r mln are least dimension and least 
radius of gyration respectively. 


In the case of short column c^ is taken as 
one and thus no reduction m permissible stresses is 
required. In the present work equation (3.93) Las 


been used. 
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3.8.4 REINPORCEMEITS : 

In general, the longitudinal steel reinforcement 
m a concrete compression member can be m between 
0.8 per cent to 8 per cent of the gross cross-sectional 
area. However, because of practical difficulties like 
accommodation of reinforcement and to avoid excessive 
cracking the percentage of reinforcement is limited 
to 4 only. In the present work the lower and upper 
limits for longitudinal reinforcement are taken as 
0.8 per cent and 4 per cent respectively. For 
spiral reinforcement these values are taken as 0.3 per 
cent and 2.3 per cent of the gross cross-section. 

3.8.5 PERMISSIBLE LOADS : 

For a column with helical reinforcement the actual 
axial load is not to exceed P which is given by 

P = Q 1 A k + c 2 A sc + 2 °3 "H) (3.95) 

In addition to this, the sum of the terms Ajr 
and 2 c .^ A K is also nor to exceed 0.5 F„ A„. 

0 D C C 

3.8.6 COLUMN SUBJECT TO COMBINED AXIAL AND BENDING FORCES 

In order to design a column on no crack basis, 
the column section and its reinforcement are proportioned 
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such, that the compressive stresses satisfy the inter- 


action formula 

OcTe + 'W^cb i 1 


(3.96) 


where , 


/■ 


= P/ f A n + m (A gc + 2 A h ) 


and 


) cb 


m/z 


where, Z is equivalent modulus of the column cross- 
section and is given as 


- M 

Z = 1 + (n-1)A 30 (D ool -2 C Cov ) 2 /4,' 2 


J col/64 


D 


col. 


= 11 D ool/32 + V ( D oor 2 °ooy)^ 2 W 


Furthermore, the resultant tensile stress m concrete 

is not to exceed 35 per cent and 25 per cent of the 

resultant compressive stress for biaxial and uniaxial 
bending respectively. In the present work the column is 
subjected to uniaxial bending only. 

. / / * 

<? tb = M/Z < 0.25 n>' ob +h c ) +4 0 


(3.97) 
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3.8,7 COLUMU SUBJECT TO DIRECT TENSION: 


If column is subject to occasional tensile 
force ) its diameter and reinforcements are 

designed to safely carry the same as follows: 
Considering the column section to be uncracked, 
the resistive force (T c0 2 _ u ) offered by the section 
is given by: 


'col u 


A col + ( m “ 1) ^ A sc + 2 A bV 


(3.98) 


Considering that concrete m the column cross' 
section is cracked, the permissible tension (T co 2 c ) 
is given by 


T col 0 = 1 -53 U sc f a + 2 A b f sh ) (3.99) 

The minimum of these two is taken as the 
permissible load (T co p p) f° r the column m Tension 
and this should be atleast equal to the maximum 
applied tensile force (C lQ ^). 

In the present work, the column diameter and 
reinforcements are assumed and stresses are checked so 
that it may safely with- stand the applied loads coming 
on to the columns under various conditions of loading 


as discussed earlier. 



CHAPTER 17 


FORMULATION OP THE OPTIMAL DESION PROBLEM 

4.1 INTRODUCTION : 

The optimal design problem has been formulated 
as a mathematical programming problem. A mathematical 
programming problem consists of minimizing or maximizing 
an objective function f(d), subject to certain constraints, 
g 3 (d) < , = , or > b^; 3 = 1 , 2, m, where 

d is a n-dimensional vector of design variables. 


4.2 LESION VARIABLES: 

In the present work, the design variables 
considered are as follows: 

i) Internal diameter of the silo (D), 

ii) Thickness of the vertical wall (T ) , 
m) Thickness of the hopper bottom (T^), 

iv) Height of the rmg-girder (H ) , 

v) Width of the ring-girder (B ), 

vi ) Diameter of supporting columns (D col ) , 

and reinforcement as percentage of cross-section as: 
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vn) 

Vlll) 

ix) 

x) 

XI ) 

XII) 

XIII) 

xiv) 

XV ) 


d 


Hoop reinforcement for vertical wall (A ^ ) » 

Vertical reinforcement for vertical wall C-A- Sw 2 ^’ 

0 / \ 

Meridional reinforcement for hopper bottom (A ) > 

Hoop reinforcement for hopper bottom (Ag^)? 

Upper longitudinal reinforcement for ring-girder 

(Agr-j ) » 

lower longitudinal reinforcement for ring-girder 

U sr2 )> 

Stirrups or shear reinforcement for ring-girder 

^ A sr3^ ’ 

Spiral or helical reinforcement for column (A_ rt _), 

SC O 

longitudinal reinforcement for column (A gc -j_). 

So the design vector becomes: 



/ 





! 

I 


T w 

% 

H r 
B r 
j D col 
: A S wi 
A sw2 
' t A sh1 
j A sh2 
A sr1 
A sr2 
A sr' 7 * 


A, 

A 


SOS 

SCI 


/ 


I 


( 4 . 1 ) 


and 



TO 


4.3 OBJECTIVE EJECTION: 

The objective or merit function provide a 
measure of usefulness or the basis of compansion 
between the acceptable designs. In general for the 
civil engineering structures, cost is of paramount impor- 
tance. Therefore, the cost of the silo super structure 
has been considered as the objective function. In order 
to account for the lift, the cost of materials for every 
metre lift above a height H~ c from the base of the 
column has been considered to be an additional 1 per cent 
more than the cost of materials upto the height H^ c . 

The total cost of the silo super structure m terms of 
the design variable can be expressed as follows: 

(Cost)rj, = (Cost)^ + (Cost) R + (Cost)g (4.2) 


where, (Cost)^ 
(Cost ) R 

(Cost ) s 


represents the cost of supporting columns, 
represents the cost of ring-girder, and 

represents the cost of silo super structure 
which consists of the roof, the vertical 
wall and the hopper bottom. 


C c i use * hsi ; Hf0 + (H c - H f0 ) (1 + r lcc )} 


4 - fy ) 

SOS ~3Cl', } 


or; (Cost ) m 


h - C r1 



/ J 


+ O r 


( 1+ p lcr ) 11 & + T ,J H B /l + C (A +a + a $ 

lcr w x x \ rl ' A sr1 +A sr2 +A sr3j- 


+ °o Vf (l+P i cruf) (^en)T^ ?*»* 


Sine. 


(1 + C. 2 A ) l. 

nl smm' 

+ t (1 + h=v> (D + T w> *w H„(1 + hfhwsl * A.) A 


+ J (1 + p . ) ii°ES£ 


H 


ich ; 2 T h Sina ^ + C r1 ^sh2 + A sh2^* 


(4.3) 


where 


H roof 

= Tan 0^ (p-p 

r ropen 

H hop 

= Tan a (D-D )/2 

K op en J/ 

^smm 

= 0.003 

°r1 

= Costal 


)/2 


P. 

ICC 


= (H 0 - H fc> he/2 


P. 


lcr (H c " H fc + H r/2^ P ic 


'icruf (H c + H r + H vw “ H fc + H rootf/2> P ic 



P 

1C7 

= »0 

+ H- 
r 

H fc + H vw/2^ P ic 

and 

p 

ich 

■ < H c 

+ H r 

" H fc “ H hop/ 2 ) P ic 

where , 

ho 

= 0.01 




and Cost 


ratio of cost for steel and concrete 
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4.4 CONSTRAINTS: 

A design is considered to be safe only when 
it performs within certain limits. This imposes certain 
constraints on the behaviour parameters of the structure 
which m turn are the functions of the design variables. 

Such constraints are termed as behaviour constraints. In 
the present problem there are 24 constraints related 
with the bounds on the various kinds of stresses developing 
on the structure. These are grouped as 1 to 24 m the 
following listing. Over and above the behaviour 
constraints, there are adways limitations appearing on 
the variables due to constructional and other practical 
considerations. These are termed as side constraints. 

In the present problem, all the fifteen design variables 
have been bounded between the lower and the upper limits 
thus, contributing to a total of 30 side constraints. 

These are grouped as 25 to 54 m the following listing. 

In order to check whether the silo or deep bin condition 
exist or not, a limit on H/D ratio is checked which 
correspond to the 55th constraint m the list. Moreover, 
to ensure the stability of the silo wall, I.S. Specifi- 
cations^^ put a limit on the wall thickness related with 
the diameter of the silo. This has been incorporated as the 
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56th constraint, 'To facilitate to use of design 
table for the torsional properties of the ring-girder 
cross-section (refer section 5 t 2.p, the upper and lower 
limit on the ratio of the two sides of the ring— girder 
has been specified. These form the 57th and 58th 
constraints. Finally, to ensure the sufficiency of 
reinforcements m the ring-girder as upper longitudinal, 
lower longitudinal and stirrups the last three constraints 
(59 to 61) are incorporated as discussed m section 
3.7.4. The -various constraint equaxions and expressions 
arranged m order are as follws. 


1 . 


Tensile stress m concrete of the vertical wall 
(without shrinkage consideration) due to hoop 
force is given as 



1 


/ 


P, 

des 

100 T w ( 1 +(m-u A sw1 ) 



(4.4) 


2, Tensile' stress m concreie of the vertical wall 

(with shrinkage consideration) due to hoop force 
is given as 


P des +6 s E s A sw 1 100 P w < 1 . 333 /^ 
100 T W (1 + (m-1) A gw1 ) ~ Ut 


(4.5) 
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3) Tensile stress in hoop steel subject to hoop force 
(considering concrete to be cracked) is given as 


' des 


si 


A * 100 T 
swl w 


-si 


(4.6) 


4) 


The net tensile stress in concrete of the vertical 
wall m case of empty silo with lateral load is 
given as 


S sh - S sw < ^tb 


(4.7) 


5) The total compressive stress m concrete of the 

vertical wall m case of silo full without lateral 
load is given as 


S sw + °sg 


< 0.15 f. 


(4.8) 


6) The total compressive stress m concrete of the 
vertical wall m case of silo full with lateral 
load is given as 


S + S + S . < 0.2 f 

sw sg sh — c 


(4.9) 


7) Tensile stress m concrete of the hopper bottom 
subject to meridinel force (without shrinkage 
consideration is given as 



FM 


des 


100 T h ( 1 + (m- 1 ) A sh1 ) 



( 4 . 10 ) 


/ 
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8 ) 


9) 


Tensile stress in concrete of the hopper bottom subject 
to mendinal force (with shrinkage consideration) 
is given as 


u t 


M des + e s E s A sh1 100 \ 
100 I h (1+(m-1) A sh1 ) 


< 1.353 0 t 


(4.11 ) 


o 

Tensile stress m mendinal steel subject to 

o , , 

meridinal force (considering concrete to be cracked) 


is given as 


FU 


des 


si 


< f. 


A sh1 100 


( 4 . 12 ) 


L h 


10) Tensile stress m concrete of the hopper bottom 

subject to hoop force (without shrinkage consideration) 


is given as 


£ t 


FH 


des 


< < 


100 T h ( 1+(m-1 ) A gh2 ) 


(4.13) 


11) Tensile stress m concrete of the hopper bottom 

subject to hoop force with shrinkage consideration 
is given as 
, / 




FH des + e s E s A sh2 100 


100 


L h 


( 1 + (m-1) A sh2 ) 


.<1.333^.(4.14) 



76 


12) Tensile stress m hoop steel of the hopper bottom 
subject to hoop force (considering concrete to be 
cracked) is given as 


PH 


des 


•si 


A 


sh2 


100 T 


< f 


si 


h 


(4.15) 


13) The sum of the shear stresses m the ring-girder 
due to sheer end torsion at any section is given as 

q + q' < 4 <^ s (4.16) 

* 4 

14) The compressive stress c (stresses and 0 ) 

m concrete of the ring-girder subject to combined 
bending and thrust at mid-span points with large 
(small) eccentricity is (are) calculated and the 
following relation is checked. 


c <_ 
cb 



m case with large eccentricity v 

| ( 4 . 1 7 ) 

< 1 in case with small eccentricity jj 


15) The net tensile stress t(S^.) m steel (concrete) 
of the ring-girder subject to combined bending 
and thrust at mid- span point with large (small) 
eccentricity is calculated and the following 
relation is checked. 
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t _< f £ ln case with large eccentricity 

/ * > s' 

- v c < 6 th in CPSe fCLth small eccentricity 

(4.18) 

^ / 

16) The compressive stress c (stresses (j and. (j 0 ) 

m concrete of the ring-giider subject to combined 
bending end thrust at support points with large 
(small) eccentricity is (are) calculated and the 
following relation is checked 


c < 




+ 



m cose with large eccentricity 
< 1 m case with small eccentricity 


f(4.19) 

5 


17) The not tensile stress t(S^) m steel (concrete) of 
the rmg-girder subject to combined bending and 
thrust at support points with large (small) 
eccentricity is calculated end the following 
relation is checked 


t < f 


s2 



m case with large 

- tC < ibb in oase 


eccentricity 


with small eccentricity 


5 

5(4.20) 


18) 


For no tension failure of the column cross-section, 

C, and T n as discnbed m section 3.8.7 

lodt colp 

should satisfy the relation 


C lodt - T colp 


(4.21) 
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At top section the colunn is subject to a 
d' wnward load, a bending moment and a horizontal force. 
The stresses are celculoced on no era cl c "oasis and the 
following relations are checked. 


19) 

C 1 A k 

+ 2 )> 

0.5 

O 

o 

(4.22) 

20) 

^ + 
(Tc 

^ < ! 

C cb 



(4.23) 

21) 


J!kr * 

-<? G < 0.25 ( 

^cb 

+ c o) 

(4.24) 


At bottom section of the column also the 
stresses are calculated on no crack basis corresponding 
to the forces acting at that section and the following 
relations are checked. 


22) 

C 1 ^k + 2 * °3 ^b ^ 3.5 ? c A c 

(4.25) 

23) 

t j 

^c^^c + ^cb^^c"b — ^ 

(4.26) 

24) 

tb “ - 0,2:3 ^ ^ 3b 

(4.27) 

25) 

"bo 

Limits on design variables: The 

following 

54) 

constraints specify the lower and 

upper bound on 


the various design variables, 


(Wi (Ui’ 1=1 - 2 15 (4 - 28) 
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where, (<i mixl )-, and (dj. ) 1 are the lower and upper 
limits on the ith design variable and are given in Table 4.1. 

TABLE 4.1 LIMITS Oil DESIGN VARIABLES 


Variable 

Unit s 

jO -er Jj.in.it; 

Upper limit 

L 

Hi 

1 .0 

9.0 

T w 

m 

0.15 

0.3 

T h 

Hi 

0.15 

0.3 

H r 

HL 

0.3 

1.8 

B 

r 

m 

0.3 

1.8 

D col 

,n 

0.4 

1.2 

A sw1 

per cent of 
cross-section 

0.3 

2.3 

^“sw2 

- do - 

0.3 

2.3 

^shl 

- do - 

0.3 

2.3 

A sh2 

- do - 

0.3 

2.3 

^srl 

- do - 

0.5 

2.3 

^sr2 

- do - 

0.3 

2.3 

^sr3 

- do - 

0.3 

2.3 

^scs 

- do - 

0.3 

2.3 

A scl 

- do - 

0.8 

4.0 
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55) Deep silo criterion is given by the equation 

H/D > 1.5 (4.29) 

56) Limitation on tie wall thickness, for a silo 
more than 6mm diamexer, is given as 

T w > 15 + (D-6 )/l . 2 (4.30) 

57) Limitations on the relative dimensions of the 
to 

58) rmg-girder are given as 

1.0 < Ratio < 4.0 (4.31) 

where , 

Maximum of H and B 

Ratio = - — — — 

Minimum of E r and B r 

59) The total reinforcement (percentage of cross-section) 
to 

61 ) required for the rmg-girder as upper and lower 

longitudinal reinforcement ond stirrups as 
discussed m section 3.7.4 are calculated as A sr ^ r > 
il sr2r and A sr5r . The relations 


^srlr 

< 

^srl 

(4.32) 

^ fc sr2r 

< 

A sr2 

(4.33) 

■^sr3r 

/* 

\ 

A sr3 

(4.34) 


are checked for sufficiency of reinforcements. 
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4.5 NORMALIZATION OF DESIGN VARIABLES: 

All the design variables are normalized such, 
that they lie in a range of 0 to 1 . Thus if d , 

(^inirPi alld ^max^i d - eno ' tos actual value, the 
lo^er and the upper bound for tne ith design variable, then 
the normalized variable (D n ) is given as 




a 




(Span) . 


i=1 ,2, . . . , 15 


(4.35) 


where , (Span) 1 = (d^^)^ - 


Such a normalization puts equal weight age to all design 
variables during the course of search for minimum design. 
The normalized variables are non-dinensional. 


4.6 NORMALIZATION OF CONSTRAINT EQUATIONS : 

The constraint equations are also normalised 
with an identical at™ as tuvt for design variables. 
Typrcal examples of normalizing a behaviour and a 
side constraint is given below: 

A behaviour constraint is typically represented 



as 


(4.36) 
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which can be rewritten as 



> 0 


dividing both the sides nyQ + - e get the normalized 

constraint as 


*> 


= 1 



> 0 


(4.37) 


A side constraint is +ypically represented as 

1 — ~ ^^rnax^l (4*38) 

Taking first part of the inequality and dividing by 
(Span)^ we have 

~ ^ q 

( Sp an ) vj 

or, (3 n ) 1 > 0 (4.39) 

Considering other part of the inequality and dividing 
dv (Span) ^ we get 

(Span ) 1 


> 0 



S3 


The above equation can De rewritten as : 


^ (3 ma:Pl “ 

(Span) ^ 



or, 1 - (h n ) 1 > 0 


(4.4-0) 


Both equations (4.39) and (4.40) are satisfied, 
if we satisfy only one equation given by 

V ( Vi L 1 - > 0 

In this way all the constraints subject to upper 
and lower bounds are reduced to only one. There are 
16 sich constraint equations m the present problem. 

The normalization of the constraints as 
discussed above effectively has reduced the number of const- 
raint equations to 45. 

4.7 NATURE OF THE OPTICAL DESIGN PROBLEM: 

Since the objective function and the behaviour 
constraints are non-linear functions of the design 
variables, the optimal design problem of the silo as 
formulated above is a non-linear programming problem. 

The method for seeking the solution of this non-linear 
programming problem is discussed m the next chapter. 



CHAPTER V 


OPTIMUM SEEKING- METHOE 

5.1 INTRODUCTION: 

The optimum design of the silo structure has been 
formulated in the previous chapter as a non— linear mathe- 
matical programming problem which is of the form: 

Emd d, such that, f(d) is a minimum, 
subject to g^(d) >0, j - 1,2, . . . . m. 

If n be the order of the design vector, then m the 
n-dimensional search space, there will m general be a 
number of design points (set of design variables) which 
will satisfy all the constraint equations and this chapter 
deals with the algorithm of picking up the best design 
corresponding to which the value of f (d) is minimum. 

Many algorithms are available to seek the solution 
of a non-linear mathematical programming problem. These 
can be classified as : 

i) Direct methods m which the constraint equations 
are handled explicitly. Well-known algorithms of this 
class are 1) Zoutendijk’s Method, and 2) Rosen’s Gradient 
Projection Method. 
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ii) Indirect methods m which the constraint equations 
are implicitly satisfied. L class of algorithms called 
Penalty Function Method is well-developed m this category. 
This has been widely used with success on a large class of 
problems. 

Picking a method is by and large an art. This mainly 
depends on the type of the problem and the tools readily 
available m hand. For the present work, the Interior 
Penalty Function Method which has advantage of giving all 
points m the feasible or acceptable region, and thus 
advantageous from engineering point of view, has been 
used. The SUMT algorithm which pertains to this class 
has been described m short. 


5.2 SEQUENTIAL UNCONSTRAINED MINIMIZATION TECHNIQUE 
The technique of optimal solution seeking, 

/ 1 O 'j 

developed by Fiacco and McCormick' , known as the 
Sequential Unconstrained Minimization Technique (SUM) is 
an Interior Penalty Function Method. In this method the 
constrained problem is transformed into an unconstrained 
problem by appending the constraints to the objective 
function through a penalty parameter r. This gives some 
sort of a high imaginary wall at the constraint surface 
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which does not permit crossing-over into the infeasible 
domain. The resulting function, F (d, r), known as penalty 
function is given as 

- m 1 

P(a, r) = f(a) + r ^ — "Tt") (5 - 0 

This function F(d, r), is minimized as an unconstrained 
— > 

function of d, for a fixed value of r and the value of 
r is reduced sequentially. The sequence of minima so 
obtained converge to the constrained minimum of the problem 
as r -*■ 0. 

5.3 UNCONSTRAINED MINIMIZATION 

For unconstrained minimization, a number of well 
developed algorithms ^ .g., Powell’s Method, Davidon- 
Fletcher-Powell Method etc.) exist m the literature and 
are documented m standard text books on mathematical 
programming. Fox^"^ has very rightly mentioned that 
picking a method is largely an art, although, to be sure, 
the better one's knowledge of the theory, the better he 
is at the art. 

In the present work, the Davidon-Fletcher-Powell 
(DFP) method has been used for unconstrained minimization 
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and is discussed m the following section . 


5.3.1 David on-Fletcher-Powell Method: 

The Davidon-Fletcher-Powell (DFP) method ^ < “* 5 
also known as Variable Metric Method is a very powerful 
and well-used method. The algorithm for this method can 
be summarised as follows : 



d + a* 

1 q. 



(5.2) 


where, d^ and d + .| are the design vectors at the qth 
and (q+1 )th iterations, 

is the move vector, and 

is the minimizing step length lh. direction 
computed by linear minimization. 


The move vector at any point 


d^ is obtained as 


where , 
and 





is a positive definite symmetric matrix 
=V5’(d q ) 


(3.3) 


is the gradient vector at 



(5.4) 
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is updated after every iteration using the formula 


H 


q+1 


= H + a* 
1 1 




(HQ) (HQ) 

v Q. .V v q 


T 



(5.5) 


m which Q^ = G-g + ^ 



(5.6) 


to start with is taken to be the identity matrix. 

5.4 LINEAR MINIMIZATION 

Eor ill-behaved functions, of the type of penalty 

function which emerges in the minimum cost design of R.C. 
( 22 ) 

Silo, Golden Section Search Technique v ' is supposed 
to be the best. This method requires only function 
evaluations. The interval within which the minima is 
bracketed is successively reduced by *618034 m this 
method.. In one linear minimization ten internal reductions 
are performed which requires eleven function evaluations. 
The algorithm runs as follows : 

i 

Let do be the starting point and the linear 
minimization is to be performed along direction, S, 
along which the minima is bracketed within an interval 
(Dist ) . |Hie steps are : 
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1 ) Initialize 
Take 


K1 = 0, K2 = 0 and ¥ = 0, 
R = 0.618034 
d = do, and 


Denote Po = P(do) 


2) Set N = M + 1 

= R times Dist 

a g = Dist - 

(Lj = d + S 
and &2 = d + ccg S 

3) If equals 1, go to step 4j otherwise continue 
Denote P^ = P(d^ ) . 

go 

4) If K 2 e 9 u als 1, A to step 5, otherwise continue. 
Denote P 2 = P (dg). 

5) If N equals 10, go to step 8, otherwise continue 

If P^ is greater than go to step 6, 

otherwise continue. 

Set cT = d 2 

V p, 

K .= 0 

1 

K 2 = go to step 7 



yu 

6) Se-t F.j = P 2 


K 2 = 0 


7) 

Set 

Dist = a.| , go to step 2 

8) 

Set 

P m in = Minimum of F^ 

and F 2 


If 

P mm e U^als F 1 , go to 

step 9, otherwise continue 


a* 

= (42 ± - ao i )/S 1 


9) 

a* 

= (41 ^ - do i )/S 1 



The a* is the desired minimizing step length. 

5.5 GRADIENT EVALUATION 

As the design-analysis cycle for the present silo 
problem uses some table, it is not possible to write expre- 
ssions for gradient of the function F explicitly m terms 
of design variables. Finite difference scheme has been 
used to calculate the gradient. Backward difference and 
central difference schemes were tried and m both cases, 
the gradient values were approximately equal. However, 
m the present work, central difference scheme has been 
used. 



For penalty function 


F(a, r) = f(d) + r ? =77 

0=1 S 0 (4) 


Pi^ck Sd 1 = minimum of d /lOOO and 0,0001, and set 


m =1 



and )d2t = 



Denot 


e f., = f (ai ), f 2 = f(d2) 


g1 D = Sj (dl) and g2^ = (d2) 


0 = 1 , 2 , 


m g2 — gl 

Then, G l( J = / g ^ /(2 S^) 

0=1 ( s 0 ) 

( 5 t 7 ) 


where, G- 1 is tne ith component of (d), 

i = 1 , 2, .... n. 



CHAPTER VI 


ILLUSTRATIVE EXAMPLE 


6.1 DESIGN DATA 


A 400 t capacity circular reinforced concrete 
silo with conical hopper bottom supported on six columns 
and used for the storage of wheat has been considered. 

The optimum design obtained through the computer programme 
developed m the present work, based on earlier discussions 
has been compared with the field design which has also 
been performed on the basis of I.S. Specifications and 
has been recently executed at Ashok Nagar (Madhya Pradesh - 
India). The data considered for the problem are as follows 

The angle which hopper makes with the horizontal (a )=45 
' Diameter of outlet opening (-D 0 p en ) = 30 cm 

Clearance upto silo opening (V Q -p) = 3.15 m 

Effective length factor for column = 0.80 


Since, the stored material is wheat, the following values 
have been used m numerical computations as suggested by 



I.S. Code 



*2 

Unit weight of the sot red material (w) = 850 kg/mr 

Angle of internal friction for stored material 

($) = 28° 

Angle of wall friction m filling (<J)* ) =*75<|) 

Angle of wall friction m emptying (<j)' ) = 

Pressure ratio m filling =0.5 

Pressure ratio m emptying (*X e ) =1.0 

Cost of concrete per cubic metre (C ) = Rs. 150/- 
Cost ratio for steel and concrete (cost ) = 75 

The additional cost of concrete for every metre 
increase m the constructional height above the 
the height for the fixed cost (P ) = 1 per cent 

Height for fixed cost (H fc ) = 3 m 

The IBii 704-4-1401 computer at IIT Kanpur has been 
used. The results obtained have been presented in tabular 
and graphical form. The numerical experience and the 
behaviour of the problem have been discussed. 


6.2 RESULTS AND DISCUSSIONS: 

The optimum design of 400 t capacity R.C. Silo 
has been obtained first by starting from six different 
starting point and choosing: 
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is followed with, the constant reduction factor m r. 

This of cour se is at the expense of computer time. While 
the typical time taken (exclusive of compilation time) 
to arrive at the optimal design with varying reduction 
factor is around 4.25 minutes, it is only 3-.1 minutes 
with constant reduction factor. Henceforth, only varying 
reduction factor to r as indicated above is used. 

The results shown m Table 6.1 indicate the existence 
of local minima m the problem. Therefore, a set of 
another nine varymgly different starting points were 
considered and optimal designs were obtained. The 
results so obtained are given m Table 6.2. These 
results further established the existence of local minima 
m the problem. 1 plot of all the fifteen optimum 
design points obtained with varying reduction factor 
m the penalty parameter versus the optimum cost of 
silo structure is shown m Fig. 13. This plot is tending 
to be asymptotic to a total minimum cost of Rs. 27,000/ 

A comparisionof the design variables corresponding to 
optimum design points having a total minimum cost of 
about Rs. 27,000/- does not show that the design is 
almost similar and from this it can be said with limited 
confidence that global minimum has not yet been obtained. 

At the optimum points, it is observed that the cost cf 
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supporting structure i.e., column and rxng-girder 
varies from 32 per cent to 41 per cent of the total 
cost while that of the super structure varies between 
59 per cent to 68 per cent. The height to diameter ratio 
(H/D) varies between 2.7 to 5.0. 

It is observed that constraint equations no. 37 

and 38 are bounded at all the optimal designs. These 

represent respectively the lower bound on the hoop and 

vertical reinforcement m the vertical wall of the silo. 

The thickness of the vertical wall also tends to near 

minimum which is reflected from constraint equation no. 32. 

Prom this it can be concluded that the optimum design 

tends towards a silo diameter which results into near 

minimum wall thickness of the vertical wall with minimum 

reinforcement both m hoop and vertical directions. This 

is justified from the fact that with the decrease m the 

silo diameter, the pressure (both m the lateral and 

vertical direction) decreases at any point along the 

depth of the stored material (the variation of lateral 

pressure with depth for different diameters has been 

plotted m Pig. 14). With the decreased lateral pressure, 

the hoop force further decreases resulting into reduced 

as 

thickness of wa.il/yrell as reduced hoop reinforcement. 
Furthermore, for a given capacity the decrease in 
diameter requires increase m height, and larger the 
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height of vertical wall more severe will be the stresses 

at the base of the wall due to lateral load, self weight 

of the wall and due to friction between the wall and 

the stored material. The optimal design, therefore, 

a 

tries to tend to/diameter which results into minimum 
vertical reinforcement with near minimum thickness to 
resist the stresses m vertical wall. As such the optimum 
design is being governed by the lower bounds imposed on 
these variables. The other variables also adjust accordingly. 

It has been observed that the cost of the silo 
structure is quite sensitive to silo-diameter and to the 
depth of ring-girder. A plot of these variables versus 
the total minimum cost for all the 15 design points 
obtained is shown m Figs. 15 and 16. These plots 
do not establish any easy mathematical correlation and 
only go to establish that the problem has a number of 
local minima. 

The cost of a 400t capacity R.C. Silo field 
design executed at Ashok Hagar is Rs. 31813/-. The 
lowest cost optimum design obtained during the present 
work is Rs. 27123/-. This leads to a reduction of 
about 15 per cent 111 total cost of the material used. 

The field design and the lowest cost optimal design are 
shown m Table 6.3. However it may not be out of place 
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to mention here that all the optimum designs obtained 
m the present work have the total cost lower than the 
cost of the field design. 

The IBM 7044-1 401 computer takes about 5 minutes 

i 

m the compilation of the programme. Additional 3 to 4 ^ 
minutes are required for one complete run (i.e., starting 
from a starting point and reaching to a optimal point) 
depending upon the choice of the starting point. 


6.3 CONCLUSIONS : 

The programme developed in the present work is 
capable of improving the design of R.C. Silo. The automated 
optimum design as formulated in the present work has a 
number of local minima. Even with a numerical experiment- 
ation of fifteen design points it has not been possible 
to reach the global minimum. This obviously requires more 
intense search. However, a trend is observed that the 
optimum design tends towards a silo dimeter which 
results into near minimum wall thickness and the minimum 
hoop and vertical reinforcements m the silo wall as well 
as near minimum hopper thickness and reinforcements m 
it. 
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TABLE 6.1 OPTIMUM LESISIT WITH VARIOUS STAKPIflO POINTS 
(Showing Infltssttsfeof Reduction m r also) 




Parti- 

culars 

units 

Starting ^5 aLMSfe 

point CRF m T/RF in 

(D 

Start- 

ing 

point 

(2) _ 

Optimal Point 

CRP in VHP 
r in 

r 

Start- 

ing 

Point 

(3) 

Optimal Point 

CHP In VBF 
r in 

r 



1 

2 

3 

1 

2 

5 

1 

2 

3 

h 

m 

5»00 

4.942 

5.151 

4.700 

4.749 

5.056 

5.200 

■iX ■■ m im ■■ 

4.857 

5.042 

T 

w 

m 

0.165 

0,164 

0*166 

0.165 

0.159 

0.165 

0.165 

0.155 

0.161 

\ 

m 

0.165 

0,157 

0,163 

0.165 

0,164 

0*160 

0.165 

0.155 

0.162 

H r 

m 

1.300 

1.258 

1.261 

1.000 

1.046 

0.964 

1*500 

um 

1,491 

B r 

m 

0.800 

O.64O 

0.538 

0.800 

0.673 

0.634 

0,800 

0,559 

0.502 

B col 

m 

0.700 

0.597 

0.603 

0.700 

0.603 

0.610 

0,700 

0.596 

0.591 

A sw1 1 

3 . of 
c/s 

0.550 

0.309 

0.301 

0.550 

0.307 

0.306 

0.550 

0.301 

0.303 

A sw2 

t 9 

0.550 

0.310 

0.301 

0.550 

0.314 

0.310 

0.550 

0,302 

0.301 

A sh1 

It 

0.550 

0.498 

0.408 

0.550 

0.506 

0.485 

0.550 

0.397 

0.450 

A sh2 

It 

0.550 

0.431 

0.560 

0.550 

0.465 

0.490 

0.550 

0.368 

0.558 

A sr1 

It 

1.000 

1.262 

1.085 

1.000 

1.058 

0.996 

1,000 

1.205 

0,890 

A sr2 

It 

1.000 

1.046 

1.057 

1.000 

1.083 

1.104 

1,000 

1.071 

1.059 

A sr3 

It 

1.900 

1*807 

1.758 

1.900 

2.033 

1.918 

1.900 

1.976 

2.030 

A scs 

II 

0.550 

0.725 

0.452 

0.550 

0.509 

0.558 

0.550 

0.439 

0.527 

A sol 

It 

2.200 

2.205 

2.252 

2.200 

2.528 

2.525 

2.200 

2.346 

2.303 

H c 

m 

4.200 

4.213 

4.314 

4.350 

4.329 

4.564 

4.100 

3.955 

4.030 

H 

vw 

m 

23.798 

24.364 

22.407 

26.966 

26.403 25.269 

21.983 

25.234 23.403 

e/d 

m/ in 

4.760 

4.930 

4.350 

5.737 

5.559 

4.602 

4.228 

5.195 

4.642 

Tcost)“ 

Rs* 

55964 ~ 

28986 

271 23 

34774 “ 

“ 28224 " 

”27216“ 

I7042" 

27737 

“27396 

(Cost) g 

P.of C T 

58,200 

60.8 

62.900 

66.0 

63.1 

63.8 

55.9 

60.2 

61.0 

(Co8t) supp " 

41. 800 

39.2 

37.1 

34.0 

36.9 

36.2 

41.1 

39.8 

39.0 

Bounded 

constraints 


30,32,33 

37,38,39 

40 

30,32, 33 ; 
37,38,39 
44 

9 

30,32,33 

37,38,59 

40,44 

30,32,35 

57,58,39 

40 

30,52, 

33r38, 

59,40, 

44 

30,32, 

33,37, 

38,39 


Continued Table 6.1 
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Parti- ni _ 
cuiar tg 

For St. Pt- 

4 

For 

St Ft 

5 


For St. Pt. & 


i 

2 

3 

1 

2 


3 

1 

2 

3 

D m 

5.0 

4.942 

5.102 

4.7 

4.729 


4. 908 

5.300 

4.998 

5.195 

T . “ 

0.165 

0*1 66 

0.167 

0.165 

0.159 


0.161 

0.165 

0.163 

0.170 

\ m 

0.165 

0.161 

0,155 

0.165 

0.158 


0.158 

0.165 

0.154 

0.163 

H r m 

1.500 

1 ,461 

1.515 

1.500 

1.477 


1.509 

1.500 

1.481 

1.490 

B m 

r 

0.800 

0.590 

0.494 

0.800 

0.566 


0.552 

0,800 

0.537 

0.547 

D . m 

col 

0.700 

0.596 

0.592 

0.700 

0.594 


0.596 

0.700 

0.594 

0.587 

A sw1 P */ f 
c/s 

0.550 

0.307 

0.301 

0.550 

0.310 


0.307 

0.550 

0.322 

0.304 

A " 

W 

0.550 

0.304 

0.310 

0.550 

0.300 


0.336 

0.550 

0.304 

0.311 

A " 

A sh1 

0.550 

0.453 

0.511 

0.550 

0.420 


0.460 

0.550 

0.403 

0.515 

A 11 

sh2 

0.550 

0.408 

0.532 

0.550 

0.434 


0.441 

0.550 

0.624 

0.525 

A » 

A sr1 

1.000 

0.830 

0.758 

1.000 

0.884 


0.999 

1.000 

1.235 

I.140 

A " 

sr2 

1,000 

1.234 

1.108 

1.000 

1.062 


1,051 

1.000 

0.891 

0.850 

A " 

a sx 3 

1.900 

1.998 

1.964 

1.900 

2.095 


1.823 

1.900 

1.969 

1.879 

A " 

scs 

0.550 

0.517 

0.722 

0.550 

0.621 


0.492 

0.550 

0.826 

0.603 

A " 

sol 

2.200 

2.247 

2.117 

2,200 

2,342 


2.258 

2.200 

2.020 

2.349 

H m 

c 

4.000 

4.011 

4.036 

3.850 

3.887 


3.945 

4.150 

4.01 8 

4.107 

H m 

vw 

00 

On 

r- 

• 

CM 

24.363 

22.846 

26.966 26.638 


24.712 

21,152 

23.814 

22.023 

H/D m/m 

4.760 

4.929 

4.478 

5.737 

5.633 


5.035 

3.991 

4.764 

4.239 

(Cost)^ Rs. 

37281 

29204 

27720 

37846 

28700 


28387 

36957 

28435 

28717 

(Cost) B of q. 57.7 

S 

60.5 

67.6 

60.5 

6 2.0 


62.0 

55.0 

61.0 

60.7 

(COSt)gjjpp „ 

42.3 

39.5 

37.4 

39.5 

38.0 


38.0 

45,0 

39.0 

39.3 

Eoanded 


30,32, 

33 30,33, 

37 

30,32, 

33 

30,32, 

33 

30,32,33 

30,33,37 

Constraints 

xs — — — — ™ 


37,38, 

40,44 

39 37,38, 

39 

37,38,39 

40,43 

37,38, 

40,44 

39 

37,38,39 

38,39,40 
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TABLE 6.2 OPTIMUM DESIGN MTH V1RI0US STARTING POIUTS 


Parti- 

TJ 

Starting Opt. 

St. 

Opt. 

St. 

Opt, 

St. 

Opt. 

culars 

n 

Point 

Pt. 

Pt. 

Pt. 

Pt. 

Pt. 

Pt. 

Pt. 


X 

ts 

i 


2 


3 


4 


D 

in 

5.500 

5.315 

4.500 

4.959 

5.300 

5.228 

4.000 

5.156 

T 

w 

m 

0.165 

0,171 

0.165 

0.165 

0.175 

0,170 

0.165 

0.173 

T h 

m 

0.165 

0,160 

0.165 

0.157 

0.190 

0.189 

0.165 

0.156 

H 

r 

m 

1.500 

1.469 

1,200 

1.177 

1.000 

1.025 

1.500 

1.41 8 

B 

r 

m 

0.800 

0.465 

0.800 

0.579 

0.750 

0.527 

0.800 

0,496 

S col 

m 

0.700 

0.597 

0.700 

0.604 

0.600 

0.630 

0.700 

0.597 

A sw1 

P.of 

0.550 

0.309 

0.550 

O.314 

0.600 

0.316 

0.550 

0.355 

c/s 









A sw2 

si 

0.550 

0.335 

0.550 

0.306 

0.580 

0.333 

0.550 

0.518 

A sh1 

it 

0.550 

0.447 

0.550 

0.500 

0.600 

0.521 

0.550 

0.564 

A sli2 


0.550 

0.350 

0.550 

0.405 

0.600 

0.525 

0.550 

0,551 

A sr1 

n 

1.000 

1.195 

1.000 

0,981 

1.400 

1.444 

1.000 

0.789 

A sr2 

it 

1,000 

0.745 

1.000 

1.163 

1.800 

1.861 

1.000 

1.003 

A sr3 

TJ 

1.900 

2,023 

1.900 

1.823 

2.100 

2.039 

1.900 

1.908 

A scs 

t! 

0.550 

0.488 

0,550 

0.480 

0.600 

0.588 

0.550 

0,566 

*scl 

n 

2.200 

2.250 

2.200 

2.323 

1.900 

1,982 

2.2 

2.214 

H 

c 

in 

4.250 

4.189 

4.050 

4.303 

4.650 

4.589 

3.500 

4.150 

H 

vw 

m 

19.622 

21.032 

29.437 

24.193 

21.152 

21.745 

37.514 

22.542 

h/d 

m/m 

5.568 

3,957 

6.542 

4.878 

3.991 

4.159 

9.328 

4.389 

Total 

Cost 

Cm 

Rs, 

36852 

27506 

56636 

27606 

37440 

28481 

40501 

29906 

T 

Sup. 

Cost 

F. of 

C T 

53.4 

62.7 

66.2 

63.6 

59.8 

62.4 

67.8 

67.0 

Supp. 










Str. 

n 

46.6 

37.3 

33.8 

36.4 

40.2 

37.6 

32.2 

33.0 

Cost 










Bounded 


30,33,35 


30,32,33 


30,37,38 


30,33,37,38 

Constraints 


37,38,39 


37,38,39 


39,40 




40,44 40,44 


Continued next page 



Contd. Table 6.2 
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Parti- 

culars 

u 

n 

i 

ts 

s t .pt. 

5 

Opt.Pt. 

St.Pt. 

6 

Opt.Pt 

D 

in 

6.00 

5.688 

7.00 

6.049 

T 

w 

m 

0.180 

0.184 

0.200 

0.187 

T h 

m 

0.220 

0.219 

0.230 

0.226 

H r 

m 

1.400 

1.440 

1.500 

1.545 

B r 

m 

0.700 

0.408 

0.700 

0*386 

B col 

m 

0.750 

0,652 

0.800 

0.666 

A sw1 

F. of 
c/s 

0.750 

0.482 

0.800 

0.328 

A sw2 

If 

0.6 00 

0.336 

0.600 

0.315 

A sh1 

tf 

0.750 

0.714 

0.750 

0.639 

' ri sh2 

If 

1.300 

1.258 

1.400 

1.280 

A sr1 

If 

1.400 

1.366 

1.500 

1.493 

A sr2 

n 

1.500 

1.481 

1.500 

1.450 

A sr3 

IT 

1.400 

1.414 

1.300 

1.484 

A scs 

IT 

0.800 

0.755 

0.600 

0.583 

A scl 

II 

2.00 

2.111 

2.500 

2.679 

H c 

m 

4.600 

4.404 

5.000 

4*479 

H w 

m 

16.441 

18.328 

11.992 

16.170 

h/d 

m/m 

2.7^0 

3.222 

1.713 

2.673 

Total 
Cost C, 

Rs* 

T 

40371 

30733 

46648 

31006 

Sup. 

Cost 

P.of 

C T 

55.7 

64.I 

51.1 

59.6 

Supp.Str." 

Cost 

44.3 

35.9 

48.9 

40.4 

Bounded 

Constraints 


15,19,30 

35,37,38 

15,19,35 

37,38 

15,19,33 

37,38 


St.Pt. 

7 

Opt.Pt 

St.Pt 

8 

. Opt.Pt St. Opt. 

Pt. Pfc. 

9 

6.400 

5.815 

5.700 

5.501 

6.000 5.602 

0.191 

0.185 

0.180 

0.179 

0.180 0.181 

0.210 

0.210 

0.195 

0.194 

0.195 0.194 

1-350 

1.393 

1.000 

0.985 

1.000 1.004 

0.500 

0.392 

0.850 

0.624 

0.850 0.624 

0.700 

0.649 

0,700 

0.628 

0.700 0.627 

0.600 

0.435 

0.700 

0.463 

0.700 0.472 

0.600 

0.432 

0.600 

0.363 

0.600 0.371 

0.580 

0.576 

0.600 

0.577 

0.600 0.577 

1.000 

0.992 

0.600 

0.581 

0.600 0.581 

1.500 

1 .466 

1.500 

1.461 

1.500 1.463 

1.500 

1.455 

2.000 

1.957 

2.000 1.958 

1.900 

1.846 

1.800 

1.764 

1.800 1.766 

0.800 

0.762 

0.800 

0.763 

0.800 0.763 

1.800 

1.960 

1.600 

1.931 

1.600 1.946 

4.850 

4.515 

4.850 

4.765 

5.00 4.797 

14.412 

17.524 

18.250 19.613 16.441 18. 906 

2.252 

3.014 

3.202 

3.565 

2.740 3.375 

36506 

31176 

39095 

30497 38946 31857 

59*4 

64,0 

57.6 

61.5 

55.6 60.7 

40,6 

36.0 

42.4 

38.5 

44.4 39-3 


15,27,30 

30,37,38 30.37 


35,37,38 


38 
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TABLE 6.3 M DESIGN JLKD OPTIMUM DESIGF 


Particulars 

Units 

Field Design* 

Minm. of opt. 




Design** 

D 

m 

7.000 

5.151 

T w 

m 

0.175 

0.166 

% 

m 

0.250 

0.163 

H r 

m 

1.500 

1.261 

B r 

m 

0.500 

0.538 

D col 

m 

0.5642 

0.603 

A sw1 

P. Of c/s 

0.865 

0.301 

A sw2 

1 f 

0.301 

0.301 

A sh1 

I I 

0.805 

0.408 

A sh2 

f f 

1.424 

0.560 

A sr1 

1 t 

0.400 

1.085 

A s*2 

t f 

0.600 

1.057 

A sr3 

f ! 

1.000 

1.785 

A scs 

f l 

0.463 

0.452 

A scl 

f r 

3.700 

2.252 

H c 

m 

5.000 

4.314 

L 

V w 

m 

11.992 

22.407 

H/D 

m/m 

1.713 

4.350 

Total Oost 

Crp Rs * 

31813 

i ! 

1 ! 

1 

1 1 

I 

! oj 

1 

I 

Sup , Go st 

P. of Cy 

65 

62.9 

Supp . Str. 

P. of Cm 

35 

37.1 

Co st 

Bounded Constraints 

18,19.50,38, 
41 j 45 

30,32,33,37,38, 
39, 44 

* Modified 

field design 

** Giving : 

minimum cost out of all the 

15 optimum designs 
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6.4 SCOPE FOR FUTURE WORK: 

The silo structure can be considered as a 
shell of revolution and rigorous analysis of the shell 
under the loads can be embedded m the programme of 
automated optimum design. 

The formulation of ihe problem given m the 
present work can be slightly modified by specifying 
all those design variables which are tending to near 
minimum values as pre-det ermine d design parameters at 
these minimum values (e.g.,the reinforcements m silo 
wall and hopper bottom). This will reduce the siae of 
the problem m the design space and thus may improve 
the nature of the function to enable us to reach the 
global minimum. 
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FIG 10 GIRDER CROSS-SECTION UNDER ECCENTRIC LOADING 




Silo - Capacity = 4 00 1 
Starting Point No 2 



FIG 12 INFLUENCE OF PENALTY PARAMETER OVER COST 
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FIG 15 TOTAL MINIMUM COST VERSUS SILO-DIAMETER 



FIG 16 TOTAL MINIMUM COST VERSUS GIRDER - DEPT! 
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APPENDIX 


COMPUTER PROGRAMME 


The programme for the present work was written for 
IBM 7044 Computer in Portran IV language. The complete 
programme consists of: 

(A) a main programme named INPUT which reads the input 
data, generates the normalized variables, goes to 
unconstrained minimization routine, reduces the 
value of penalty parameter, checks the convergence and 
prints the output, and 

(B) the other eight subroutines namely: 

I ) DEPOWL which provides the general working 
of the DPP method and is called from INPUT, 

II ) GRAD which calculates the gradient as per central 
difference scheme and is called from DFPOWD, 

in) HASH which updates the Hassian Matrix and is 
called from DPP01D, 

iv) GOLDEN which performs the linear minimization 
as per Golden Section Search Technique and is 
called from DEPOWL, 

iv) PENAL which converts the normalized variables 
to the actual variables, goes to the analysis- 

design subroutine, calculates the constraints 

related with upper and lower limits on the 
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design variables and finally gives the 
penalty function value , 

vi ) ANLSIS which calculates the cost of the silo 
structure and all the constraints other than 
those related with upper and lower bounds on design 
variables and is called from PENAL, 

vn) SPECTR which gives the average acceleration as 
per I.S. Specif ic?rion''"and is called from the 
ANLSIS , and 

vui) TORS ON which gives the torsional properties for 

a rectangular cross-section as per Table 3.2 and is 
also called from ANLSIS. 

The complete listing of these sub-routines are 


as follows: 





C-Giuu FORTRAN SOURCE LIST 

ISN SOURCE STATEMENT 

0 SI8FTC INPUT 

CC ' . 

C AUTOMATED OPTIMAL DESIGN OF 1 

C REINFORCED CONCRETE 

C SILO . j 

r ■ j 

1 DIMENSION C{ 5),DN(15),DMIN(15),SPAN(i5),GD(45) 

2 DATA : i , MC , NON/ i 5,45 ,3 / 

3 DATA jPSit GPS2/.-.4 :Ci» -j.O. >011001/ 

4 DATA DNIN/-U. 0,2*0. 15, 2*0.3,0.4,8*0.003,0.008/ 

5 DATA SPAN/3. ; ,2*. 15 ,2*1.5, 0.3,8*0.02,0.032 / 

6 DATA GCONCM, $ IGC8, S IGC, $ IGTB, S'IGT, $IG-S»S IGBON , S I G 8R / 
i E ; • f . ' , 5 - • -* , j / *t, 1 2 * 0 , IT «v> ,8.0 ,40 .0 / 

7 0 AT A F S , F S H , F 3 1 , F S 2 /1300.0, LCOO.0, 1000.0, 1250.0 / 

10 DATA c S / ,1 . iO + 0 6 / 

11 DATA W,FI , WCtOOPcN, TTA / 850.0, 28. 0,2400.0* 0.3 ,45 .0 . / 

.2 DATA WSEPF » WINDPR, E S£?F , SI SC03, FZ3 , W BETA/O. 7, 1 50. » 1 . 5 , .04, .3, 1.0/ 

V *3 DATA 8C0V , CCOV, ELF ACT , VC LEAR/ 0.04, C* 04,0. 8, 3.15/ 

14 DATA WTQP, TMT,C RAT I C , TTARUF, THR , ROPEN/lOO* , .0 , 75. ,5 . , .08 , • 3 / 

15 DATA CC / ,5 .3. / 

* i6 P I = 4 . 0 * A T A N (1.0) 

17 CALL FL UN (.0000) 

ID A AMU C R= 3 -3C 0 . /GC ONCM 

2 i I F C AMODR.GT. 18. i ) AMCDR-iS .0 

24 A M 0 0 R I = A M 0 D R - 1 . 0 

25 cC = 55/ AKGOR 

a 6 FI=FI*PT/.Sv.y : 

27 T T A = TT A * P 1 / 1 IU « . > 

30 TT ARUF=TTA RUF*P 1 / 13 C. 0 

31 AMU=TAN (FI) 

•32 F I F = ) . 7 5 * F 1 

33 FIE® . .6 * F I 

34 AM1JF=T AN ( F I F ) 

35 AMUS=T AN ( F IE ) 

36 ALL A F=, j, .5 

37 ALC A c= i. . J 

• 40 COMMON /NM8ER/ NFE, KGctNGS ,NMA 

■ 41 COMMON /PRG8/NC ,NCN 

42 COMMON /NORML/ DM IN ,SPAN 

43 COMMON /CONCR/GCONC M, SIGCB , SI GC, S IGTB,SIGT,SIGS,SIG80N,SIGBR, A MO DR) 

44 COMMON /STEEL/ FS, F SH , FSi , FS2 , ES , EC, AM0DR1 , PI 

45 COMMON /COST/FON6W , CC » CR» CAP, WTOP , TMT, W, WC , FI , DOPEN, TTA, AMUF , AMUE 

*ALDAF, ALDA E, COSTS, CCSTC,CQSTR, CRAT 10, TTARUF,THR, ROPEN I 

hb COMMON/DTFT/BCQV ,CCCV , 6LFACT, VCLEAR, WSEP F, W INDPR , ESEPF, S I SC03, FZ3 : 

* , WBE TA 

47 COMMON /SHAPE/ H,HC 

50 111 CONTINUE : 

51 PRINT H 

52 CALL TIME(M) 

, . . ' :"■* S3 - NFE*0 DM" 

54 READ2, CAP, (D( IT, 1=1 ,N ) 

• -61 PRINT!, CAP, (D(I ), 1 = 1, N) 

66 1 FORMAT ( * OPT. COST DESIGN OF FOOD GRAIN STORAGE SILO*, 5(3H **) , /IX, 

*42 (1H~ ),///* CAPACITY®*, F10. 5,* TONNES*,//* VARIABLES ARE . .*,/* 

* MIDI A. AND (2) THICKNESS OF SILO CYLINDER, (3 ) THICKNESS OF SILO 


d14- - 



160 FORTRAN SOURCE LIST INPUT 

1SN SOURCE .STATEMENT 

♦HOPP-R, (4) HEIGHT AMD (5) WIDTH OF RING GIRDER*,/* (6) DIA.OF SUPP: 
PORTING COL. AND PER . RE INF . AS (7) HOGP AND (3) VERTICAL REINF+FGR V.: 

* WALL*,/* (9) MERID INAL AND UO) CIRCULAR REINF.FOR HOPPR. ( 11 ) UP . ' 
*AMD (ID LOWER LONG. AND (13) STIRP REINF.FOR R.G .*,/* (14) SPIRAL : 

* AND (.3) LONG. REINF.FOR COL. . *, / i X, 3F10. 5 , /IX, 3FiO. 5) : 

67 2 FURR AT { •'* Fi 0 . 5 ) 

T, 3 FORMAT ( * COM SIR AN IS ARE* , / 10 (6F16.7,/) ) ‘ 

71 4 FORMAT ( * STARTING PCI NT IN INFEASI RLE DOMAIN*) : 

72 3 FORMAT ( /* COST T=*,E13.7,* R=*,E15.7,* FUNC.V.=*» E15.7,*C0ST OF SU: 

*P . 6 , CO L » R . 6 . * , 3 ..13. £ ) 

73 6 FORMAT ( /5 EX, *NUM8ER OF F UNC . H VAL . * » 1 5 , /5uX» * NUMBER OF GRAD . EV AL .* , . 

*15,/3EX,*MUMLCR OF GOLDEN SER* , 1 5 , /50X , *NUMB ER OF MATRIX ASS*, 15/3.' 

74 7 FORMAT ( * COST IMP. IN LAST I TER . *, £ 15 . 7 , / * PENAL PART* ,E15 . 7 , / 

** TOTAL , SUP.S. , SUP. C. , R.G. COSTS ARE*,4E13.5 ) 

75 8 FORMAT ( * >:RR?=* , 5 15 .7,* . GOING FOR NEXT CY. WITH R=*,£-15.7) I 

76 9 FORMAT ( * ERR .-.=*, cl 5 .7 , * . GOING FOR NEXT CY.WITH R=*,E15.7) 

77 10 FORMAT ( * AT THE ABOVE POINT THF RATIO OF HT.TO 01 A . =* , F8 .3 , * VERT 

*. WALL F ;; I G HT = * , F 8 . 3 ,* COL. HEIGHT**, FB. 3) 

1-EA ii FORMAT { •. Hi) 

t >.'•* DO i.5 I-i,N 

102 15 ON ( I ) = ( U { I ) - OMIN ( I ) )/ SPAN { I ) 

104 CALL PEHAL(F0,FP,FT,N,DN,1.0,GD) 

105 PR I M3 , ( GO ( I ) , 1 = 1, NC) 

U2 HuR = 1 .5 / ( 1 . -GO ( 1 ) ) 

113 PR I Ml 0 , HDR, H, HC 

114 COST E=FO 

115 1 F ( FT..LT . 1 . • 6+25 ) GOTO i? 

."120 ‘A - : PRINTS : 

12X ERR 1=1.0 

122 E RR 2=1. A 

123 GOTO 4 ■ 

124 17 ' DO 18 1=1, NC 

125 TFCGD( I) . LT.Q.1 ) FP= FP~1 . O/GD ( I ) +10.0 

13C 18 CONTINUEE ■ A - . . A . A'A/A; a :: A , "A ' 

132 NPR = AL‘.jG10(F0/FP) 

133 . R =5.0*10. 0**NPR 

134 RCHECK = 'R/1.0E + 07 

135 NCH5KR= .> 

136 FT=FO+R*FP 

137 PRIMS, FO, R* FT, COSTS, COSTC ,COSTR 

140 20 CALL DFPOWLl N,DN,R, FT) 

141 CALL PENAL ( FOjFP, FT ,N, DN, R, GD) 

142 PRINT3, C GD ( I ) ,1=1, NC) 

, ■ 147 " AaA'-AEA PR 1 NTS , F0 , R, FT, COSTS, COSTC, COSTR - 

150 PRINTS, MFE, NGE, NGS, NMA 

151 NCH£KR=HCHEKR+1 

15 2 GOTO ( 21 , 21 ,21 ,22 , 22 ,22, 22 *23 ,23, 23 ,23, 23,23,23, 23 ,23 » 23 , 23 ) , NCHEKR 

153 21 R=R/2. 0 

154 GQTC 2+ 

155 22 R=R/5. 0 

156 GOTC 24 

157 23 R=R/10. 0 

168 .24 CONTINUE 

161 DO 25 1=1, N . 

162 25 0 ( I ) =D MI N ( I ) +DN ( I ) * SPAN ( I ) 



mgi6; 


FORTRAN SOURCE LIST INPUT 


I so 


SOURCE STATEMENT 

1 64 


PRINTS} (0(1) ,I=i,N) 

% -> 

4. * 4. 


HDR = ..« 3/ (1 .O-GO ( 1 5 ) 

IT 2 


PRINT! , i HDR»Hj HC 

173 


1RR1=( COST, ,-FONtW) / FONSW 

L 1 4 


CUbT FOHfr.W 

17 5 


1RR2= ( FT-FOfLiW) / FUN EW 

176 


IF (SRR1.LT.EPS 2.) GOT C3 

2?/- 1 


PR INTS ? - :: RR2}R 

202 


GOTC 35 

203 

3 0 

I Ft ERR. .LT.tPSl )G0TC4 

26 6 


PR IN, T9, t. RR.,} R 


3 5 

CONTINUE 

210 


I .Ft R.Lf.RCHECK) G0T040 


C 

I F { R.LO. i » c-..7)G0TC4. 

213 


PR INTI 98 ■ 

2 i A 


GOTC 21 ■' 

Zi 5 

4 0 

PRINT?, !RR„, 1RR2, COST., COSTS, COSTC ,COSTR 

2x6 


PRINT! 9 7 » H DR 

2;. 7 


PR I NT;. 9 9 

2 2 U 

1 97 

FORMAT ( /* AT THE DPT. COST POIN THE RATIO 

221 

1 98 

FORMAT ( l X,*2U(1H-) , /) 

222 

199 

FORMA T(//1X, 40 ( 3H##— ) ,// ) 

223 


GOTC 111 

224 

2 00 

STOP 

22 5 


2ND 


C5G160 


IBMAP ASSEMBLY INPUT 


FOR ABOVE ASSEMBLY 


FORTRAN SOURCE LIST 


C cb i 6 J 

ISN SOURCE STATEMENT 

0 Si BFTC DFPOWL 

1 SUBROUTINE DFPOWL { H ,D , R, FO ) 

C 

C D A V I 0 0 M - F L £ T C H c R-PUWELL METHOD 

C OPTIMUM SEEKING ALGORITHM 

C 

2 DIMENSION h{ .6) , S { 1. 6 > ,G(i6) , Gl (16) ,Q( 16) ,HJ (16,16) 

3 D I HENS I ON GO ( 43 ) , DM IN ( 1 5 ) , SP A N ( i 5 ) 

4 COMMON /NMB'dR/ NFS , NGE , NGS, NMA 

5 COMMON /PROB/NCtHCM 

6 COMMON /NGRML/ DM1 N, SPAN 

7 COMMON /COoCR/GCONCM, SI GCB ♦ S I GC , SI GTB , SIGT, SI GS , S IGBON, SI GBR , AMuDfl 

IN COMMON /STEEL/ FS, F SB, FS1 , FS2 EES , EC, AMODRl, PI 1 

11 COMMON /COST/FONEW , CC ,CR » CAP , WTOP, TMT , W» WC , FI , DQPEN , TT A , AMUF , AMIUE ■ 
*ALDAF, ALDA E , COSTS * CCS TC » COST RjCR AT 10 f TTARUFf THR »ROPEN 

12 COMMON/r.vTFT/BCOV»CCCV , EL F ACT , VCLEA R» WSEPF » WINDPR , ESEPF , 5 1 SC03 , FZ3 

* » WBt 7 A 

13 NFt- 

, . 14 1 G c — 0 

15 0GS=0 

* 6 N M A = 

17 1 FORMAT!* FUNC . VALUE =* » £13 .7 , /* AT STARTING VECTOR*, /2 ( 8 £15 . 7 , / ) ) 

20 2 FORMAT (* GRADIENT*, /5 (8216. 7,/ ) ) 

21 3 FORMAT!* DIRECTION*, /5{8E16. 7, /) ) ! 

22 4 FORMAT ( * WITH ST£P=*,£X5.7,* FUNC. VALUE AND COSTS ARE*,4E15.7) 

23 5 FORMAT { * FUNC .VALUE =* ,£15.7,* WITH STEP=* ,E15.7 , * AT NEW POINT *,/ 

*2 { 8 S 15 .7 , / ) ) 

24 7 FORMAT ( * FORCED CON V.*, 13,* ITERATIONS PERFORMED* 

*/* STEP CONV. CHECKED LI AND L3 ARE *,215,* RESPECTIVELY*) 

25 8 FORMAT (* NORMAL CONV.*, 15,* ITERATIONS*) 

26 9 FORMAT ( * X . . . FU KC .ERR =*,£16.7,* GOING TO RESTART ... X*) 

27 I T ER=1 

30 . . L 1 = C ' ■ ■ ■ 

31 L3=0 

32 CALL PENAL (FOl, FP» FTi ,N, D, R, GD) 

33 IF{ FTl .GT.L.Ofc 25IRETURN 

36 Fu=FT! 

37 PRINTljFO, (DU) ,1 = 1 ,N) 

44 CALL GRAD ( D, N, G » R ) 

C INITIALIZATION OF H J MATRIX, DIR. S AND PRESERVATION OF G AS Gl 

45 21 DO 30 1=1, N 

46 00 25 J= 1 , N 

47 25 H J ( I , J ) =0 . 0 

51 30 H J { I , I ) = 1 . * 

53 31 DO 41 I = 1 , N 

54 Gl ( I )=G ( I ) 

55 S( 11=0.0 

56 DO 40 J=1,N 

57 40 S ( I J =S ( I l-HJ ( I , J ) *G { J ) 

62 FQLC-FQ 

63 SS=O.D 

64 SSG=0.0 

65 DO 45 I = I , N 

66 SS=SS+SU)*S(1) 

67 45 SSG=SSG+S ( I )*G( I ) 



FORTRAN SOURCE LIST DFPOWL 

ISO SOURCE ' STATEMENT 

71 $$- SQR T ( SS ) 

C CHOCK AS S TRASPGSF.OcL F/2.F 

72 £RR=ABS£ SSG/2./FO) 

73 I F { ERR .LT.i. OS- 1 0 )G CTGSO 

76 IFISSC.LT. , / 5 G 0 TO 5 C 

C OTHERWISE DIRECTION IS NOT USEFUL 

i'. i GU »C 2 . 

102 50 CONTI M U 

C DIRECTION IS USEFUL 

C INTERVAL FOR SEARCH IS BEING FIXED 

C NO CAR I ABLE IS ASSUMED TO TAKE NEGATIVE VALUE 

103 01 ST =2 • 0 

- A- JU 5 ■; l -it \ 

105 S(I)=SU)/3S 

106 1 F ( S ( I ) . GL * 0 .0) GOTO 52 

ill LIST :M=-D( I )/SU ) 

1*2 I F ( C 1 S TEN . LT .01 ST ) D IS T = DI STEM 

115 52 CONTINUE • 

C as PER NEED DIST MAY BE RtDUSDED BY DIVD SUITABLY BY 10 ORIOO 

C IN SIR. OPT. PROS. SPECIALLY 

117 01 ST=D I S T / IS . 

120 I F ( A BS { DIST) . GT . i . 0 E-3 ) G0TQ5 5 

123 ST 3P=0 . J 

124 L2 = L3+i .. 

1 25 GOTO So 

126 55 CALL F 130 { G, N , S » DI S T, STEP , FQ , R ) 

127 ; I F t STEP. LT .1 . 05-35 ) L3=L3+I 

132 I F ( STEP • LT . 1 *u2-G5 ) LI— Ll+1 

135 56 CONTINUE : 

CC 

136 IFtLl.GE.S -IGQTQ.70 

141 1 F £ L3. GE • 2 )G0T07:» 

144 I F ( S TE P . L E . 0 • 1 E-3 5 ) G0T02I 

CC 

C ITERATION LIMIT MAY BE REDUCED FOR SAVING TIME 

147 I T R = 1 Tt R + i 

150 IF{ ITER.GT.2 ?)G0TG7C 

CC 

153 CALL GRAD(D,N,G,R) 

154 12=0 

C THIS CHECKS CONV.OF D..SIGN VECTOR 

C STEP LIMIT MAY BE CHANGED AS PER NEED 

CC 

155 I F { A BS ( $ T E P ) . L £ . 1 .0 E- 35 } GO TO 6 5 
CC 

160 00 57 1=1, N 

161 IF LABS {STEP* ST I ) ) . L T. I .0E-05 ) G0T057 

164 12=12+1 

165 57 CONTINUE 

167 IF(L2.6Q.0)G0T075 

172 65 DO 67 1=1, N 

173 67 Q( 13 =G ( I ) — G ! £ I J 

175 CALL HASH ( N, HJ, S , Q, STEP ) 

176 GOTO 31 • 

177 70 PRINT7,ITER,L1,L3 


CcGI 60 


FORTRAN SOURCE LIST DFPOWL 

I Si i 


SOURCE STATEMENT 



tRR = A8S( ( FC- FOL D ) / F C ) 

2 u 1 


GOT C 85 . 

20 2 

75 

cRR=ABS( { FO— FOLD) /FC ) 

203 


IFCERR.LT. 0.0001) GQ TO 80 

C ' . 6 


PRINT9, ERR 

20 7 


GOTO 21 

22 0 

80 

PRINTS , ITER 

2 i, i. 

85 

PRINTS, FO, STEP, (DID, 1=1, N3 

216 


RETURN 

227 


END 

C t G 1 6 0 


IBMAP ASSEMBLY DFPOWL 

NO MESSAGES 

FOR above ASSEMBLY 



ucbitK; FORTRAN SOURCE LIST 

ISM SOURCE STATEMENT 

0 $1 BFTC GRAD 

1 SUBROUTINE GRAD < D,N ,G,R) 

C 

c GRADIENT EVALUATION BY CENTRAL DIFFERENCE 

: c 

2 DIMENSION D( 16) ,GU<£) 

3 DIMENSION GO ( 45 ) » DM IN M5 ) » SPAN ( 1 5 ) >GD1 (45 ) » GD2 ( 45 ) 

4 COMMON /NMBER/ NFE » NGE , NGS , NMA 

5 COMMON /PRQB/NC , NCN 

6 COMMON /NORM!/ DMIN » SPAN 

7 COMMON / CONCR/GCONC M» SIGCB t S IGC » SI GTB » SIGT * SIGS » S IGBON » SI GBR » A MOD 

10 COMMON /STEEL/ FS»FSH T FS1»FS2»ES»EC»AM0DR1,PI 

11 COMMON /COST /FONEW» CC ,CR , CAP , WTOP > TMT » W t WC » F I » DO PEN , TTA, AMUF, AMUE 
*ALDAF, ALDAS , COSTS » C CSTC ,COSTR, CRAT 10 1 TTARUF, THR , ROPEN 

12 COMMON/DTFT/BCOVfCCCVj EL FACT i VCLEAR, WSEPF , WINDPR, ESEPF, S I SC03 , FZ3 
* j W D E T A . 

13 NGc=NGE+I 

14 CALL PENAL(FO,FP,FT»N,D,R,GD) 

15 DO IT 1=1, N 

16 ST = AMIN1 (O.OQOl »D( I )/ 1000,0) 

17 D ( I ) =D ( I ) -ST 

20 CALL PtNAL(FQI,FPI, FT_ , N , D , R , GDI ) 

C MIDDLE DIFFERENCE 

21 D( I )=D(I J+ST+ST 

22 CALL PENAL(F02» FPZ t FT2»N,D,R,GD2) 

23 G { I ) = ( F02-F01 ) / ( 2.0 *ST ) 

24 DO 5 J=l,NC 

25 5 G( I )=G( I)-R*( (GD2( J J-GOM J) )/(2.0*ST) )/(GD(I)**2) 

27 0 ( I ) =D ( I ) -ST 

30 10 CONTINUE 

32 RETURN 

33 END ' 

CEG16Q IBHAP ASSEMBLY GRAD 


NO MESSAGES FOR ABOVE ASSEMBLY 


C.cGl 60 


ISN 

0 

1 


SOURCE STATEMENT 


FORTRAN SOURCE LIST 


CEG160 


■EIBF1C HA SM 

SUBROUTINE HASM (N, H J, 5, Q, STEP ) 


2 

3 

4 

5 


C 

C 

c 

c 

c 


6 

7 

10 

12- 


13 


14 


15 


16 


17 

10 

22 


23 


24 


25 


26 

20 

M::; 


32 


33 

30 

36 


37 



DFP 


UPDATING OF THE HASIAN MATRIX FOR 

Ul MENS ION HJ { 16 , 16 ) »S { 16 ) , Q( 16 ) 

~nw wn?, 1 ? fSi HWi { 1 6 > * HNN ( 16 » 16 ) , HNN (16,16) 

UURRUN /NM6ER/ NF£ , NGt ,NGS » NMA 

NMA =NM A+l 

REASSEMBLY OF HASSIAN MATRIX 
( JQ + 1) = ( JQ}+ (iMOJ+CNC) 

OQ=G 0+1-GQ 

(MQ ) =A LFASTAR*( SQ*S£) /SQ*QQ 

{ NQ)=- { ( JQ)*QQ)*( JO }*QQ/QQ( JQ)QQ 

HMO = :: . j ' 

DO 5 I = ], f N 
HNi { 11 = 0,0 
DO i 1 = 1, N 
HMQ=HMO+$( I )*Q( I ) 

DO 10 J=I,N 

HNI ( I )=HNi ( I )+HJ ( J, I)*Q( j) 

HMN(I, J]=S( t)*s( J) 

CONTINUE* ' 

HND = 0« 0 

DO 20 I=I,N 

HND = HND+HN1( I )*.Q-{ I ) 

DO 20 1 . J = 1 , N 

HNNU, J)*HNl(n*HNJ, (J) 

DO 30 1=1, N 
DO 30 J=1,N 

ReURN i=HJ(I,J) +STEP * HMN{ I»J>/HMD-HNN(I, JT/HND 

END 

IBMAP ASSEMBLY HASM 


NO MESSAGES FOR ABOVE ASSEMBLY 



FORTRAN SOURCE LIST 


C tbi 6U 

ISN SOURCE STATEMENT 


0 fcIBFTC GOLDEN 

1 SUBROUTINE FI 8.0 t D* N »S» 01 ST * STEP* FO *R ) 

c 

c gulden section search tech. for linear minmization 

c 

2 DIMENSION DU6 ) , S ( 1 6 ) , Di ( 16 ) , D2 ( 16 ) , DO ( 16 ) 

2 DI MENS ION GO ( 45 ) , DM IN { 15 ) * SPAN { 15 ) 

4 COMMON /NMBER/ NFE , AGE ,NGS, NMA 

5 COMMON /PROB/NC i NCN ? 

6 COMMON /NORML/ DM IN »SPAN : f 

7 COMMON /CONCR/GCQNCMt SIGC8 , $ I GC , SI GTB , SI GT, SIGS, SIG8GN ,S IG8R, Af 

10 COMMON /STEEL/ FS, FSH,F$1 , FS2, ES,EC,AM0DR1,PI I 

11 COMMON /COST/FONEW.yCC |CR,CAP tWTOP,TMT, W* WC»FI,DOPEN,TTA,AMUF,Af- 
*ALDAF, ALDAS, COSTS, CCSTC,COSTR, CRAT 10, TTARUF, THR, RQPEN 

12 CQMMUN/DTFT/eCOVtCCCV rELFACT, VCLEAR, WSEPF, WINDPR, ESEPF, SISCOS, F 

* , WBETA ■ ■ f 

13 NGS=NGS+l 

14 NCHECK= r 

15 '' . . DO 40. 1 = 1, N 

16 40 DO ( I ) =0 ( I ) 

C GOLDEN SEARCH 

20 . LOGICAL ONE, TWO 

21 RI = 0 .6 18034 

22 R2 = l. 618034 

23 41 CONTINUE 

24 ONE=. FALSE. 

25 TWO=. FALSE. 

26 I TER=0 

2.7 51 I TER*I T6R+1 



10) GOTO? 5 


56 


30 

31 
.32 

33 

34 55 
36 

41 

42 

45 

46 

47 
. -52 

■ 55: 

56 . 
57:,. 58 
6 Cl- 


ALFA1=DIST/R2 
ALFA 2- 01 ST- ALFA l 
00 55 1 = 1* N 

Olil>=D(I I+ALFAi*St I) 

■02 ( I ) = 0( I )+ALFA2*S ( I) 

I F { ONE ) G0T056 

CALL PENAL (FC,FP, FI »N»D1,R,GD) 
IF {TWO )G0T057 

CALL PENAL (FC,FP,F2 ,N,02,R,GD) 
CONTINUE 

I F { FI . LT . i ,0E 20 ) GO 1059 
I F (F2. LT . 1 . 0 £ . 2 Q ) SO JOSft 
DIST =ALFA2 
GOTO: 41 ■' : 

01 ST=ALFA1 
- 


C EG 160 


r ' FORTRAN SOURCE LIST GOLDEN 

I SM 


SOURCE STATEMENT 

LI/ 


TWQ=. FALSE . 

1.0 i 

ftj 

. dist=alfai 

102 


GOT C 51 

103 

75 

FM I N = AM IN 1 { F l,F2) 

Iu4 


' DO 76. 1=1, fv 

105 


L= I 

i06 


IF { AB5 ( S ( I ) 3 .GT .0.1 E-303GOT077 

lii 

76 

CONTINUE 

113 

77 

J = L 

114 


.IF ( FMIN.GT .FOJGQTOiCu 

1X7 


1F( FMIN.EQ.F23G0T085 

122 


ST£P={ OH 1 3-DO ( I ) )/S{ n 

123 


DO 8 .4 I = 1 , N 

124 

80 

CM 13=0111) 

126 


FG = Fi 

127 


RETURN 

130 

85 

S T E P = { D 2 ( I 3 - DO ( 1 3 )/S( 1 3 

131 


DO 90 1=1,0 

132 

90 

0(1) =02 £ I 3 

134 


FG=F? ' 

135 


RETURN 

136 

100 

CONTINUE . 

137 


DO 102 1=1, N 

140 

102 

D { 1 3 =D0 ( I 3 

142 


NCHECK =NCH£CK+1 

1 143 


IF{NCHECK.GT.2) GOTO 117 

146 


IF(FMIN.gQ.Fl 3G0T01C4 

151 


DI ST= ( 02 ( L 3- DO ( L 3 3 / S ( L 3 

152 


GQTC UQ 

’ 153 

104 

01 ST = { Dll L 3-DLK L 3 )/S(L) 

154 

110 

IF ( CIST. L£ .X.OE-8 3 GCT0117 

157 


GOTC 41. 

160 

117 

STEP=0..0 - 

161 

120 

RETURN 

162 


END 

CEG160 


I SWAP ASSEMBLY GOLDEN 

NO MESSAGES 

FOR ABOVE ASSEMBLY 



SOURCE STATEMENT 


FORTRAN SOURCE LIST 


CcGi'60 

ISN 


$1 BFTC PENAL 

SUB ROUT I NL PENAL ( FO »FP, FT »N*DN»R* GD ) 
C 

C INTERIOR PENALITY FUNCTION 

C 


2 DIMENSION GQ5) , DN{ 15 ) , DMINC 15) ,SPANU5 ) ,GD(45) 

3 COMMON /NMBER/ NFE , NGE , MGS, NMA 

4 COMMON /PROB/NCtNCN 

5 COMMON / NCRML/ DMIN ,SPAN 

6 COMMON /CONCR/GCONC Mf S IGCB » SIGC » SI GTB » S I GT , SIG S » S IGBON » SI GBR » AMO Of 

7 COMMON /STEEL/ FS » F SH , FS1 , FS2, ES » EG* AMODR1, PI 

COMMON /COST /FONEW , CC , CR » CAP , WTOP , TMT , W , WC, FI , DGPEN , TTA, AMUF, AMUE 
*ALDAF, ALDA E, COSTS, CCSTC,COSTR,CR AT IO»TTARUF,THR»ROPEN 
il. COM MON /DTF T/ BCO V, CCCV , ELFACT, VC LEAR, WS EPF, WINDPR, E5EPF, SISCOS, FZ3 

* » hBE FA 

12 COMMON /SHAPc/ H,HC 

13 NFL = NF E + X 

14 PI =4 . :)# ATAM( 1.0 3 

15 DO 10 1 = 1 , n 

16 10 0 { I 3=DMINU )+DN( I )*SPAN( I ) 

20;' CALL ANLS I S t N , D , GD, FO ) 

C INCLUDES CAL. OF LOA C, STRES S , AND CONS. OTHER THAN ON +VENESS ANDSPN 

C ALSC INCLUDES COST CAL. 

CC 

CC 

21 DO 20 1=1, N 

22 II-NCN+1 

23 20 GD( II) «0N( I Hil.O-DM I)) 

25 SUM=0 . y 

26 DO 30 I = I , NC 

27 I F ( GD{ I 3 . L T . 0 .0 3G0TC40 

32 IF ( GO ( I ) . EG. 0.© ) GD ( I)=1.0E-08 

35 3u SUM = SUH+1.<./GD{ I) 

37 FP=R*SUM 


40 

FT=FO+FP 

41 

FONE W=FO 

42 

RETURN 

43 40 

FT=1.0E 

44 

RETURN 

45 

..END 


t 

CfcGl 60 




FORTRAN SOURCE LIST 


61 SO 

ISM SOURCE STATEMENT 


. $IBFTC AML SIS 

1 SUBROUTINE AWLS IS (M ,X,6, FO) 

C 

C EVALUATION - LOADS, STRESSES, COST OF THE 

C SILO 

C AMD CONSTRAINTS OTHER THAN BOUNDS ON VARIABLES 

C 

2 DIMENSION XI 15 ) , G ( 45 ) ,CLOD ( 2 ) , CLBM ( 2 1 

3 COMMON /PROB/NC , NCN 

4 COMMON /CONCR/GCONC M , SIGCB , SI GC » SI GTB > SIGT , SIGS,S I GBON , S I GBR, AMO DR 

5 COMMON /STEEL/ FS, FSH, FS1 , FS2, ES , EC, AMODR1 , PI 

6 COMMON /CO ST/FONEW , CC ,CR, CAP , WPOT , TMT, W, WC , F I , DOPEN , TTA , AMUF, AMU£,i 

*ALDAF, ALDAE , COSTS , C CSTC , COSTR , CR AT 1 0, TTARUF , THR , ROPEN ) 

7 COMMON/UT FT/BCO V, CCCV , EL FACT , VCLEAR, WSEPF, WINDPR, ESEPF, SISCO, FZF 

* ,'WBETA ■ ■ 

ID COMMON /SHAPE/ H,HC 

11 D=X ( i ) 

12 TW = X { 2 ) 

13 TH = X C 3 ) 

14 H R = X ( 4 ) 


15 BR = X { 5 ) 

16 D C 0 L = X { 6 ) 

17 ASW1=X (7) 

20. AS W 2 “X ( 8 ) 

21 ASHl=X ( 9 ) 

22 ASH2=X £ 10) 

23 ASR1=X ( ill 

24 ASR2=X (12) 

25 ASR3=X( 13) 

26 ASCS=X ( 14 ) 



0 )/ ( P 1*0*0/ 

■■■■■ 


ASCL=X (15) 

C ICCL=P I *DC0L**4/ 64.0 
ACOL=P I *DCOL*DCOL/ 4 »0 
AR=HR*BR 

RMIY=8R4HR*#3/12.0 
. RMR={D+TW) /2.0 
Rl=RMR-BR/2.0 
■ R2®RHl+BR/2» Ov. : 

WTOP=D*WPOT 

E = { HR +T W * S I iE ( T T A ) - TH ) / 2 • 0 
C SHAPE AND PRESSURE 

HHG P = { D-DO.PEN) *TAN ( TTa)/2.0 
HCLEAR=D*T AN (FI )/4.C 

H-ICAPW * ■.-/W-PI*FHGP*(D*D+DOP£N*DOPEf 

14.0 )+H CLEAR 1 ''I'D 
44 : .■ HT=H+HHOP ■ l . D - ^ T:/ : . v ■ : :' 

. ; C : .DEEP:. SILO CRITERIA 

■ '$5 - ' ■ ■ G( I ) = I .0- 1 .s*n/H ■ : - 1 

46 G(2i=1.0-0. ; 






CcGl 60 

ISN SOURCE STATEMENT 


FORTRAN SOURCE LIST ANLSIS 


133 WT2 = NH0PG-*-PI*D* { WC* { H*TW+8R*HR ) +WTCP ) +WHOPC 

134 wT=WTi 

C EARTHQUAKE AND WIND FORCES 

C DAMPING FACTOR IS TAKEN AS 5 PER. 

135 HC= VC Li: A R +HHG P -HR 

136 HCOL=HC#£LFACT 

137 AMaSS=WT/93*. «G 

J.40 AK£ C=6 . -j * ( 3 . *EC*CI C0L/HC**3 3*100.0 

141 TNAT =2 ;*PI / ( SORT ( AKASS/AKEQ ) ) 

14?.. CALL SPECTRt TNATjAVACCL) 

143 PEAR TH-E SE PF * AM ASS* (FZF*WB£TA )*A VACCL 

144 PW IND= WSCP F*WIMDPR* (D*H+HHOP* (D+DOPEN 3/2.0) 

145 HORZF = AHAXU PEARTH, PWIND) 

C VERT .W . VERT .RE INF . FOR VERT. LOADS AND HORZ. LOADS ETC. 

i.46 TFWF=AHAXl (TFWF,TFWE) 

147 STKESW=( WTOP+H*TW*WC)/(TW*< 1 . +AMODR 1* ASW2 ) ) 

150 $TRESG=TFWF/ (TW*( 1, C+aM0DR1*ASW2 ) ) 

15 1 STRE$H=2.0*H0RZF*H/ (P1*TW*(1 .0+AMODR1*ASW2 3*D**2 ) 

152 G ( 6 3 =1 * G- ( STRESH-STRESW)/(SI GT*10000 .0 ) 

153 G( 7 3 = 1 . 0 - ( STRESG+STRESW) / { 1 500 .0*GCONCM ) 

154 GI8M.0-1 STRESH+STRESW+STRESG)/ (2 COO.O*GCONCM) 

C HOPPER DESIGN 

155 PHDES=AMAXl( FM1*DM/C, FM21 

156 FHHDES = AMAXl ( FHUP1*CM/D»FHUP2 ) 


C ME RID I NAL RING 

157 AEQ=TH* { 1 . O+AMODRi* ASH1 3 



160 , 

161 

162 

-163 

164 

165 

166 

167 

170 

171 


G ( 9 3 =1 .0-FMDES/ ( A£Q*S I GT *10000. 03 

G( 1 C )=A.a~(FMDES+3.C*eS*ASHl*TW)/( AEQ*l3333.3*$IGT ) 
G(il)=i.Q-FMDES/{FSl*A SHI *10000.0) 

CIRCULAR OR HOOP RfcINF. 

AfcQ—TH* ( 1 .0+AM0DR1* ASH2 ) 

G(i2) = i..0-FHHDS$/( AEQ*SIGT*10Q00.0 3 
GU3) = 1.0-(FHHDES+3.0*ES*A$H2*TH)/ (A£Q*13333.3*SIGT3 
G(i4)=i.0-FHHDES/( F$1*ASH2*10GGO.G) 

INTER ACTION BETWEEK COLUMN AND RING GIRDER 
: TMT*EE*FM1 . 

REAL Ki»K2,K3,K4,K5,K6,K7,K8,K9 
DAT.ft Ki,K2»K3,K4»K5 »K&,K7,K8»K9/ 

4 . .083333,. 04797 ,-.01482,. 00751, -.866,-1.0, -.9549, .00151 , .2982/ 

DATA DKK,DDKK / 0.003364,1.1153 / 

HHGG=AMAXi ( BR , HR 3 
B B S S= A M I N 1 ( BR » HR 3 
RaTIO=HHGG/BBSS 

1 if fitful CT I ON ON R I KG GIRDER 3 ! ; ■ 

. il . o-i . 

CALL TOPRO (RATIO, ALFA, BETA, GAF 

+ ;HCl 

5/CICCLI 









CLG160 FORTRAN SOURCE LIST ANLSIS 


ISN 

SOURCE STATEMENT 



211 

CQLCUi=WT/6. 



2.2 

SFRG1= Kl *W f 



4 X 

SFRG3=K2*WT 



214 

BMMAX1 =K3*'WT*RMR 

+ K5*TMTA + 

TMT*RKR 

215 

BMMAX2 = K4*WT *RMR 

+ K6*TMTA + 

TMT *RMR 

216 

BMMAX3=K7*TMTA + 

TMf*RMR 


217 

TMMAXi ='..5*TMTA 




22sl T M M A X 3 ~ K 8 * WT * RM R + K 9 *TMT A 

C i ,2 ,3 . . REFER TO SUP . » MID . » MAX TORS POINTS RESPECTIVELY 

C RING SEAM DL5 . 

221 ANCR = 1 .- / ( l . j+F S2/ ( AMODR*S IGCB ) ) 

222 A J C R = 1 . ANCR/3 .0 

223 QCR=4NCR*A JCR*$ IGCB/2.0 

2.24 TMMAX=TMMAX3 ■ 

225 IF {BHMAX3.LT .0.0) E PMAX1 =AMI N1 ( BMM AX1 , BMMAX3 ) 

23.) 1 F ( E MM AX3 . GT . 0 » 0 ) 6 P H A X 2 = A M A X I ( B M M A X 2 , B M M A X 3 ) 

233 8MMAX=AMAX1( ABS t BMM AX1 ) » BMMAX2 ) 

234 QPwCj, = SFRG_ / (BR*AJCR*( HR-BCOV ) ) + TMMAXi/ ( BETA*BBSS**3 ) 

235 0PiC3=SFRG3/ (BR*AJCR*(HR-BCOV) ) + TMMAX3 / ( BETA*B BSS**3 ) 

236 G <16 ) * i . 0-AMAX1 { QPQ Cl f QPQ03 )/ (40GCG. 0*SIGS) 

237 4SRIB=ABS ( BMMAX1/ { FS2*AJCR*( HR—BCOV) ) J/lOQOQ.Q 

240 ASR2B=AQS ( BMMAX2/ ( F 52*AJCR*{ HR-BCOV ) ) J/1000U.0 

241 ASR3 S1 = S FRG1/ ( FS2* A JCR* ( HR-BCOV ) ) / 10000.0 

242 ASR3S3 = SFRG?> / { FS2*A JCR* ( HR-BCOV ) )/ 10000. 0 

243 ASR2T = TMMAX/(0.8*FS2*(BR-1.5*BC0V)*(HR-1.5*BCQV) )/10000.0 

244 A$i2T=ASR3T*(HR-U5*BCOV)/2.0 

245 ASR1R= ( ASR 1B+AS12T ) /AR 

246 aSR2R= ( ASR2B+AS 12T ) /AR 



;i*BCHR+A 


247 

250 

251 

252 

253 

254 

255 

256 

257 ' 

260 

C 

261 

262 

265 

266 

1 CoE y 

271. . 

272 

273 

275 ■ 

276 120 

277 

300 

' A O0 

■■ ■■ ; ; 


ASR3T1 = TMM AXi/ (0.8* FS2* ( BR-1 .5*BC0V ) * ( HR-1 .5*BC0V ) ) /10OOO.0 
ASR3T3=ASR3T 

ASRST=AMAX!,{ A5R3S1 + ASR3T1 , ASR3S3 +ASR3T3 ) 

ASR3R=ASRST* ( HR+BR- 3 . Q*BCQV ) / AR 

GM7 ) = 1.0— ASRiR/ASRl 

G ( 1 8 ) = 1 .(j- ASR2R/ ASR 2 

G( 19 )=1 . 0 - ASR3R/ASR3 

RCGMP=SX1*RMR 

P=RCOMP 

BCHR=BCOV /HR 

CHECK FOR eCCENTRTCTTY 
ERA=HR*BR* ( 1 .Q+AM0DR1* l ASR1+ASR2 ) ) 

I F ( BMMAX2 . LT . 0.0 ) PR INT921 , BMMAX2 
.. ^XENsBMMAX 2:/P : . ■ A. v 
IF(EXEN.GT. (HR/6.0) JG0TO120 
DES. ON NO CRACK BASIS - SMALL 
AN=(.5+AM0DR!*( ASRI*BCHR+ASR2*(1.-BCHR)1 
ERI=(I ./12.+ ( .5-AN ) =**2+AMQDRl*(A$Rl*( AN--E 
***2) )*BR*HR**3 

G( 2C ) = 1.0+ (-P/( ERA* SI GC ) ~BMMAX2*AN*HR/ ( EF 
G ( 21 ) = 1 . 0+ ( P/ ERA-BM M AX2* ( 1 . - AN) * HR/ER I 1 / 

GOTO 123 
ACQNS 
000R1*A. 


CcGl 64 

ISN SOURCE STATEMENT 


FORTRAN SOURCE LIST ANLSIS 


SOI AN»C.5 

30 2 DO 1 Li 1=1,10 

3-3 ANNEW=AC0NS/{AN**2+BCaNS*AN+CC0N$) 

3d 4 DIFF=ABS { AN-ANN'EWl 

305 AN = { AN+ANN5W }/2 .0 . i 

306 IF C (DIFF/AN) .LT.C.0C01 )G0T0122 ! 

311 121 CONTINUE 

333 122 CSTRES=P*( EXEN/HR+. 5-BCHR ) /{ „5*AN* ( 1.-BCHR-AN/3 * ) +AM0DR1*A5R1*( 1 ,-i 

8BCHR/ AN ) # ( i.-2.*BCHR) ) / ( BR*HR*1G00 C . 0 ) 

314 G{ 20 ) = l *0— CSTRSS/SI CCS 

315 G ( 2 1 ) = 1 .0- AMQDR4CST RE S* ( 1 .-BCHR-AN ) / { AN4FS2) 

316 123 IF ( BNMAX1. GT.Q.OJPR INT922, BMMAX1 

321 EXEN-AB S ( BMMAX1/P ) 

322 I F ( EXEN . GT .{ HR/6.0 ) JG0T0124 

325 AM= ( • 5+.AM0DR1* ( ASR2 48CHR+ASR14 ( 1 .-BCHR ) ) )/(l. + AMODRl*( ASR1 + ASR2) ) 

326 ERI = U./12.+ (.5-AN) **2+AM0DRl*( ASR2*( AN-BCHR)**2+ASR1*(1.-BCHR-ANJ 
4442) ) *8R*HR*43 

32 7 G( 22 ) = i..o+ (~P/( ERA*SIGC)+BMMAX1*AN*HR/(ERI*SIGCB) )/ 10000.0 

330 G(23 1 = 1.0+ (P/ERA+BMNAX1*( i. -AN) *HR/ERI ) / ( S IGTB*100'00 .0 ) 

S3 i GOT C 126 ' ■]; 

332 124 ACONS=&.*( { AMODRl^A SR2*BCHR+AM0DR* ASRI* ( 1 .-BCHR ) ) 4 ( EXEN/HR- .5 ) + ( AM 

OOORi*ASR2*BCHR**2+APODR*ASRl*(l.-BCHR)**2J ) I 

333 6C0NS«3.*EXEN/HR-1.5 

334 CCQNS=6 . 4 ( (EXEN/HR-«5i*'(ASR2*AM0DRl+ASRl*AM0DR)+ASR2*AM0DRl*BCHR+Al : 
SSRl*AMODR*fl.-.BCHR) ) 

335 AN=0 . 5 

336 DO 125 1 = 1,10 : ; • it 


337 ANNEW=ACONS/ { AN**2+EC0NS*AN+CC0NS ) 

340 DIFF=ABS( AN-ANNEW) • / '4 - ,4- : - ■ :;,42U4... 4 : ^t4.>, 

341 AN= { AN+ANNEW ) /2 .0 

342 IF( (DIFF/AN) .LT .0. 0C01 3 GOTO 1 26 

345 125 CONTINUE 

347 126 CSTRES=P*( EXEN/HR+. 5-BCHR 1/ ( „5*AN* (1.-BCHR-AN/3. ) +AMQDR1*ASR2* (1 .■ 



BBCHR/AN)4(1.-2.4BCHR ) } / ( BR4HR*10G0G.O) 

G{ 22 )=1 .O-CSTRES/SIGCB 

G(23)=i.0-AM0DR4CSTR£S4{ l.-BCHR-AN)/ ( AN4FS2) 

CONTINUE ■■■ ■■-■■ 


SUPPORTING COL. 

RATI 0=HCOL/DCQL 
IF (RAT 10. L '.15. 0 ) GO TO 13 U 
CRC=1.5-RATI0/3Q.G 
GOTO 145 
30 C R C = „ • 0 
35 Ci=CRC*SIGC 
C2=CRC*$IGC6 


as 


|ip|p 







FORTRAN SOURCE LIST ANLSIS 


C _Gl 6 . ! : 

ISN SOURCE STATEMENT 


373 01 V* AM INI ( (ACOL£F*SIGT) , ( ASCL*FS2+ 2.Q*ASCS*F$H ) *i.333*ACQL) 

374 G { 24 ) = 1 * 0— ( COLODH-C CL0D2 ) / ( DI V*10000. G ) 

cc 

-7 5 C LOC { I ) =C0L001 +COLO CH 

376 C L B M { L ) = A 8 S ! T MT A ) . 

371 GLOC ( 2 ) =CLQD { 1 ) +CLGE5W 

400 CL6M ( 2 ) = ABS ( HC* ( HA+PORZF/6.0 ) -TMTA ) * }\ 

401 GC5) = 1.0-(SIGC*((1 .-2.0*CC0V/DC0L ) **2-ASC L ) +2 . 0*FSH*ASCS ) / ( ..5*FC ) : 

402 G(26>=i.0-3.Q*CLBM( l)*ACQL£F/{5.0*CL0D(l)*Z) 

4.3 G { 27 ) = L . 0 - {CLOD ( 1 ) / (ACOLEF*S JGC ) +C LBM M ) /( Z*SI GCB ) ) / IGQGG.G 

404 G ( 2 3 ) = i . 0 - (Cl*( (1.-2. 0*CC0V/DC0L 3**2-ASCL )+2.0*C3 *ASCS) / ( . 5*FC ): 

405 G(29 )=I.0-3.0*CLBM! 2 ) *ACQLEF/ ( 5 . 0*CL0D ( 2 > *Z ) 

406 G( 3C ) = i » 0— ( CLOD { 2 ) / !ACOLEF*Ci ) + C LBM ( 2 ) / ( Z*C2) ) / 10000-0 
CC 

407 CRl =CR ATI 0-1 .0 

410 HRGCF= (D-RQPcN)*TAN(TTARUF)/2.0 

411 DMR=DM 

C PERCT. I NCR .IN COST - PIC - DUE TO INCR.IN CONSTR.HT. ABOVE HFC 

412 HFC=3.u 

413 PIC=Q.Q1 

414 PI C C = ( HG-HFC ) *P I C/2 .0 

415 PICR = ( HC+HR/2. 3-HFCl*PIC 

416 P I C V = ( HC + HR+H/2 . 0-H FC ) *P I C .... ,j 

417 P I CH= ( HC + HR-HFC-HHO F/ 2. 0 ) *P I C 

420 PICRUF= { HC+HR+H+HRQCF/2.0-HFC )*PIC 

421 CQSTC=CC*(6.*AC0L*(HFC+!HC-HFC)*(1 .+PICC) ) * ( 1. +CR1* ! ASCL+ASCS } )) 

422 CO$TR=CC*( 1. +P I CR ) * PI* ( 2. *RMR*AR* ( l.+CRl*( ASR1+ASR2+ASR3 ) ) 1 

423 COST S=CC*PI * ( ( 1 * + P I CV }* ( D+TW ) *TW*H* ( 1 .+CR1*( ASW1+ASW2 ) 1+ U.+PI CHJ *j 
*DM*TH*HHOP/SIN( TTA ) * ( I -+CR1* ( ASH1+ ASH2 ) ) ) 

* +CC*PI*(1.+PICRUF) *DMR'*'THR*HftQOF/SIN{T.TARUF.'l*Cl»O+CRl*0*.6O6}.' 

424 FQ-COSTC+CQSTR+COST S 

425 901 FORMAT (1X,FI0.2,6£17.7) 


434 
■ 435 

436 

437 

440 

441 
442 ; 

«EG160 


FORMAT (IX, 120 (1H- )) 

FORMAT ( * MAX .H0RZ. F IL MAX. EMT.AND DES . PRS. ARE* ,3E17.7 »* COR. 
1 *» 2F10.3) 

FORMAT (IX, 6E 17.7 ) 

FORMAT!* WT.QF HOP. GRAIN AND SILO-HOP *,2E17.7 
*/* MERIOINAL AND RING FORCES *,4E17.7) ■ v 

FORMAT (* MERIOINAL AND HOOP FORCES AT MIDHOPPER HT.*,4£17. 7} 
FORMAT!* WIND FACTORA AND FORCE ARE*, 3EI7.7, . 

*/* EARTHQUAKE FACTS .AND FORCE *,6E17.7) 

FORMAT!//* LARGE ECCENTRICITY*) . 7' 

FORMAT!*.. XX MI C. B. M.=*,E15.7,* X X . . *) 

FORM AT ! * . . XX SUP. B. M.=* ,615. 7,* X X . . *) 

FORMAT!* . . • - AN - *, £15 .71 - ■. 

FORM AT ( IX ,40 (3H**- ) ) 

■ocTimw .. 
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FORTRAN SOURCE LIST 


CcGI&C 

ISM SOURCE STATEMENT 

0 SIBFTC SPECTR 

1 SUBROUTINE SPECTR ( T» ACCL ) 

C SPECTRA FOR DAMPING FACTOR = 5 +ER FIG. 15 P. VI. 1.16 N.B.C.OF I . 7 C 

2 I F ( T .Gc • 0 • 4. AND . T.LE.0.8JGQT01 

i I F { T . G c . 0 . 8 . AND . T . L E . 1 • 2 ) G0T02 

ED I F { T . 3 L » 1 • 2 • AMD . T . L E ■ 1 .6 ) G0T03 

13 I F ( T .GE. I . 6. AND .T .LE.2.0) GQT04 

i-6 I F { T .GE .2 . 0. AND .T » L E. 2. 4) G0T05 

2 1 I F ( T . G> . 2 .4. AND .T. L E. 2.8} G0T06 

24 I F ( T . G T . 2 . 8 ) GOTO 8 

27 PR I NT7 » T 

30 STOP 

31 1' ACGL=I ? 5.0-5 , .0* ( T-0.4J/0.4 

32 RETURN 

33 2 ACCL=125.Q-35.0*<T-C.3)/0.4 

34 RETURN 

35 3 ACCL = 9»‘ . 22 .0* ( T-l .2)/0.4 

' 36 RETURN 

37 4 ACCL=6fj.0-U.0*(T-l .6)/0.4 

40 RETURN 

41 5 ACCL = 5 7.0-7.L*< T-2.C)/0.4 

42 RETURN 

43 6 ACCL=50.0-3.C*( T-2.4)/Q,4 

44 RETURN' 

45 7 FORMAT { * NATURAL PERIOD = *,E17.7,* IS OUT OF RANGE*) 

46 8 AC CL =47.0 

47 RETURN 

50 END 
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Mm 


FORTRAN SOURCE LIST 


TDPRQC AN, ALFA, BETA, GAMA) 
SAF PAP 

0. AND. AM .Lc.l .5 )G0T01 
5 , AND. AM .Lc.2. j JG0T02 
:.) . AN D . A M . L E . 3 . 0 ) GO TD3 
0 .AND. AM .L t • 4 . 0 } GCT04 


AMM= ( AM-1 . . ) /•;. . 5 

A L F A = •: j . 1 4 0 + 0 . 1 5 4 4 A M P 
BETA = 0 . 204+0.13 8*AMP 
GAMA=1 .000-0. 1414AM P 
RETURN 

A M M = ( AM-1 . 5) /u. 5 
ALF A=() . 294+0 .163 *AMP 
B IT A =.0 .3 46 +0.14 74AM M 
GAMA=3 . S5 9- j .0644AM P 
RETURN 

AMM= (AM-2 . 0) 

ALFA=o .457+0. 3334AMM 
BETA=0 .493+0. 3U84AMM 
GAMA = ?) . 795-0 .042 4AM P 
RETURN 
AM M= AM-3.0 
ALFA=0 .79+ 0.3334AMM 
BETA =0.801 +0.3494AMP 
GAMA =0.7 53-0. 00 S4AMM 
RETURN 

FORMAT ( 4 RATIO OF R INC 
10 F RANGE* ) 

END 
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IBLDR 


JOB ■ 000000 


LUDING IOCS 


00000 THRU 12251 
12260 
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